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FOREWORD 


This  report  describes  the  Arithmetic  Processor  (AP)  microprograms  (or  MACROS) 
developed  for  the  passive  synthesis  portion  of  the  Digital  Acoustic  Sensor 
Simulator  (DASS)  program.  These  microprograms  are  written  for  the  Arithmetic 
Processor  of  the  AN/UYS-1  Advanced  Signal  Processor  (ASP)  built  by  IBM,  Inc. 
Approximately  half  of  these  programs  are  highly  specialized  for  the  DASS  applica¬ 
tion  and  are  of  limited  general  use.  The  others  have  more  general  application  and 
should  be  of  interest  to  other  ASP  users. 


The  work  described  herein  was  sponsored  by  Naval  Air  Systems  Command, 
NAIR-370,  under  Airtask  5495490/004D/1W11 020000. 
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CHAPTER  1 
INTRODUCTION 


This  report  describes  the  Arithmetic  Processor  (AP)  microprograms  (or  MACROS) 
developed  for  the  passive  synthesis  portion  of  the  Digital  Acoustic  Sensor  Simulator 
(DASS)  program.  These  microprograms  are  written  for  the  Arithmetic  Processor  of  the 
AN/UYS-1  Advanced  Signal  Processor  (ASP)  built  by  IBM,  Inc.  Approximately  half  of 
these  programs  are  highly  specialized  for  the  DASS  application  and  are  of  limited 
general  use.  The  others  have  more  general  application  and  should  be  of  interest  to 
other  ASP  users.  A  brief  description  of  each  MACRO  is  given  below.  A  complete 
description  of  each  MACRO  is  provided  in  its  corresponding  chapter.  Chapter  format, 
for  the  most  part,  conforms  to  a  previous  IBM  publication.1  Each  chapter  contains 
a  functional  description  of  the  algorithm  implemented  by  the  mi croprogram,  a 
mathematical  description  giving  the  exact  computations  performed  by  the  microcode, 
a  description  of  the  microcode  implementation  from  the  viewpoint  of  sequencing  and 
control,  and  tabulation  of  register  usage  and  program  statistics.  In  addition, 
each  chapter  contains  detailed  charts  for: 

a)  Program  coding  and  timing 

b)  Arithmetic  Element  Controller  (AEC)  and  Arithmetic  Element  (AE) 
register  assignments 

c)  Working  Store  (WS)  layout  of  input  and  output  data 

d)  Scaling  analysis  of  input,  intermediate  and  output  data. 

The  program  coding  charts  are  a  fixed-format  representation  of  the  sequence  of 
micro-instructions  for  each  microprogram.  Each  line  on  the  chart  represents  the 
events  occurring  during  one  100  ns  microstep.  The  register  assignment  charts 
list  the  computational  elements  which  are  stored  in  each  AE  register  during  micro¬ 
program  execution.  The  Working  Store  layouts  illustrate  the  WS  data  organization 
used  in  DASS,  although  other  equivalent  layouts  could  be  used.  The  scaling  analysis 
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charts  list  each  partial  result  of  the  computations  and  give  its  fixed-point 
scaling  in  terms  of  the  assumed  position  of  the  binary  point.  An  Fx  scaling 
indicates  that  the  partial  result  is  a  32-bit  (fullword)  number  with  the  binary 
point  x  bits  to  the  left  (x>0)  or  right  (x<0)  of  the  midpoint  of  the  word  (16  bits 
from  either  end).  An  Hx  scaling  indicates  a  16-bit  (halfword)  number  with  an 
assumed  binary  point  x  bits  to  the  left  (x>0)  or  right  (x<0)  of  the  right  end 
of  the  halfword. 

Appendix  A  contains  the  source  listings  for  each  microprogram. 

BBAMP  (BROADBAND  AMPLITUDE) 

BBAMP  weights  an  array  of  complex  spectral  components  by  a  set  of  spectral 
amplitudes  and  convolves  the  result  with  the  spectrum  of  a  sine  pulse.  This  MACRO 
is  peculiar  to  DASS  and  has  limited  application  to  other  problems. 

DSCLN  (DISCRETE  LINE  SPECTRA  GENERATION) 

DSCLN  generates  the  appropriate  spectral  components  for  discrete  line 
components  of  arbitrary  frequency,  amplitude  and  phase  and  adds  these  spectral 
components  to  the  Broadband  spectrum  generated  by  BBAMP.  This  MACRO  is  also 
peculiar  to  DASS  and  has  limited  application  to  other  problems. 

SRFFT  (SINGLE  REAL  FFT) 

SRFFT  performs  the  additional  pass  or  operation  necessary,  in  addition  to  the 
basic  (complex  to  complex)  FFT  algorithm,  to  effect  transforms  between  real  (time) 
arrays  and  single  sided  complex  (frequency)  arrays.  This  MACRO  should  be  of 
general  interest  and  applicability. 

SCLA2  (SCALE  AND  ADO  TWO  ARRAYS) 

SCLA2  provides  the  combination  of  two  arrays,  X  and  Y,  to  produce  an  output 
array  Z.  Provision  is  made  for  independent  scaling  of  each  of  the  two  input 
arrays.  This  MACRO  is  of  general  applicability. 

ASSSS  (ASYNCHRONOUS  SAMPLE,  SCALE  AND  SUM) 

ASSSS  resamples  an  input  array  at  a  rate  which  is  non-integrally  related  to 
the  origin  sample  rate.  Linear  interpolation  is  performed  to  reduce  aliasing. 
Resampled  data  is  then  scaled  and  added  to  a  second  array.  This  MACRO  is  of  some, 
but  not  extensive,  general  usefulness. 
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SCLA3  (SCALE  AND  ADD  THREE  ARRAYS) 

SCLA3  provides  the  combination  of  three  arrays  X,  Y,  and  Z  to  produce  an 
output  array  W.  Provision  is  made  for  independent  scaling  of  each  of  the  three 
input  arrays.  The  MACRO  is  of  general  applicability. 

F22S3  (FILTER  TWO  AND  SUM  THREE  ARRAYS) 

F22S3  is  intended  for  filtering  each  of  two  inputs  and  adding  the  sum  of  the 
filter  outputs  to  a  third  input.  Each  of  the  two  filters  is  a  two-pole  recursive 
filter  with  independent  parameters.  This  MACRO  may  be  of  use  in  other  applications. 

DEMON  (DEMODULATED  NOISE) 

DEMON  is  used  to  generate  two  signal  components,  each  consisting  of  a 
broadband  signal  modulated  by  a  periodic  signal.  The  MACRO  additionally  combines 
the  two  components  with  a  third  input  to  produce  a  composite  output  signal.  This 
MACRO  is  of  some,  but  not  extensive,  general  usefulness. 
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CHAPTER  2 

BROADBAND  AMPLITUDE  (BBAMP)  MACRO 

INTRODUCTION 

The  Broadband  spectra  generation  algorithm  is  described  in  detail  in  a  previous 
2 

publication.  Briefly,  the  broadband  spectra  is  described  by  an  array  A  of  860 
numbers  representing  the  spectral  amplitude  at  860  frequencies  evenly  spaced 
across  the  band  of  interest.  The  860  frequencies  are  given  by: 

fn  =  fc(n  -  f  );  n  =  1,  2,  ....  860 

where  f  is  the  DFT  cell  width.  Each  specific  spectral  array  B  is  generated  by 
the  following  algorithm: 

8«t  *T<VRi  -  WW 

BI,  =■  i-  (AlRI,  -  Ai+1RIi+1) 

BRq  -  BI„  -  0 
Ag60  must  be  0 

BR.  =  BI^  =  0  for  i  greater  than  860 

where  BRi  and  BI^  are  the  real  and  imaginary  components  of  the  ith  cell  of  the 
output  spectral  array  B,  i=0  to  1023,  and  RRi  and  RI ^  are  independent  samples  from 
a  Gaussian  random  process  with  zero  mean  and  variance  of  y. 

The  DFT  (single  sided  complex  to  real)  of  the  B  array  produces  a  2048 
sample  time  waveform  segment  with  a  sine  pulse  weighting.  Successive  segments 
generated  in  this  manner  are  overlapped  by  1024  samples  and  combined  to  form  the 
output  broadband  time  waveform.  Since  the  segments  are  uncorrelated  (all  RR's, 

RI's  are  independent)  the  overlapped  segments  add  incoherently  to  yield  a  random 
function  having  the  desired  average  spectral  behavior  and  stationary  first  order 

^AVIS,  R.  H.,  "Synthesis  of  Steady-State  Signal  Components  by  an  All-Digital 
System",  N0LTR  74-215,  Naval  Ordnance  Laboratory,  (Now  Naval  Surface  Weapons  Center) 
05  Dec  1974. 
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statistics. 


FUNCTIONAL  DESCRIPTION 

BBAMP  MACRO  consists  of  four  steps  to  obtain  the  non-zero  broadband  spectrum. 
The  first  step  is  to  generate  a  "variator"  factor  which  is  a  continuous,  piecewise 
linear  function  defined  over  the  non-zero  amplitudes.  This  function  is  produced 
by  beginning  at  each  breakpoint  with  the  value  of  the  function  and  adding  an 
increment  (signed)  for  each  successive  value  until  the  next  breakpoint  is  reached. 
The  second  step  is  to  multiply  the  AMP  element  by  the  variator  value  generated 
as  described  above  to  produce  the  amplitude  for  that  spectral  cell.  This  amplitude 
is  then  used  in  the  third  step  to  scale  the  two  components  of  the  corresponding 
RN  (Random  Number)  array  element  to  produce  the  "half  cell"  frequency  components. 

t*  h 

The  final  step  is  to  form  the  output  array  element  by  linearly  combining  the  i 

s  t 

half  cell  components  with  the  previous  ((i-1)  )  half  cell  components. 


MATHEMATICAL  DESCRIPTION 


Phase  1:  CAVi  =  CAV^  +  DAV  j 

where  CAV-  .  is  the  variator  function  value  at  the  last  cell  and  DAV.  is  the 

* "  *  th  * 

increment  value  for  the  (current)  j  linear  segment.  If  i  corresponds  to  the 

St 

j+1  breakpoint, 

CAVi  =  AV.+1 

where  A V  is  the  initial  value  for  the  j+1  segment  obtained  from  the  VAR  buffer. 


Phase  2:  VA.  =  CAV.  *  Ai 


,th 


where  A^  is  the  AMP  value  for  the  i  cell,  and  VA^  is  the  variated  amplitude. 


th 


Phase  3:  BHRi  =  VAi  *  RRi 
BHI.  =  VAi  *  R ^ 

where  RR^  and  RI.  are  the  real  and  imaginary  values,  respectively,  of  the  i 
RN  element,  and  BHR^  and  BHI.  are  the  properly  scaled  components  for  the  "half- 
cell"  frequencies. 

Phase  4:  BR,  =  BHR,  -  BHR. 


■i  -—I  —'i-1 

BI.  =  BHI.  -  BHI.J _.j 


,th 


where  BR.  and  BI^  are  the  final  output  components  of  the  i  element  of  the  BB 
array. 


IMPLEMENTATION 

BBAMP  is  implemented  with  a  ten-instruction  main  loop  which  processes  one 
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pair  of  A's  (amplitudes),  R's  (complex  random  numbers)  and  B's  (output  broadband 
spectral  cells).  An  eighteen-instruction  preamble  loads  the  pipe  such  that  two 
valid  B's  are  produced  on  each  pass  through  the  main  loop.  An  eight-instruction 
segment  of  code  contained  within  the  main  loop  is  skipped  except  when  a  new  set  of 
variator  parameters  is  required  for  a  new  linear  segment.  The  variator  array  is 
stored  as  two  32-bit  word  pairs  -  the  first  word  contains  the  initial  variator 
value  for  the  i^*1  corresponding  linear  segment  in  the  high  order  16  bits,  and  a 
count  parameter  determining  the  number  of  pairs  of  amplitudes  for  which  the  segment 
is  to  apply  in  the  low  order  16  bits.  The  second  word  contains  the  increment  value 
for  the  segment.  The  count  parameter  is  interpreted  as  the  number  of  amplitude 
pairs  -  1,  and  thus  a  count  of  zero  produces  a  segment  of  length  2.  The  first 
count  parameter  is  interpreted  by  the  preamble  as  number  of  pairs  -  2;  and  there¬ 
fore,  the  minimum  length  of  the  first  segment  is  4.  The  count  is  tested  for  zero 
and  decremented  each  pass  through  the  main  loop.  The  embedded  variator  function 
change  code  is  executed  if  the  zero  test  is  met. 

The  inside  loop  is  controlled  by  BNZ8.  An  outside  loop  controlled  by  BNZ9 
allows  arrays  of  greater  than  512  points  to  be  processed.  BBAMP  is  intended  to  be 
combined  with  DSCLN  in  a  SUPER  MACRO,  and  returns  to  the  SUPER  MACRO  Code  via  BR4 
upon  completion. 

SCALING 

See  Table  2-1 . 

ARITHMETIC  ELEMENT  CONTROLLER  (AEC)  REGISTERS 
See  Table  2-2. 


ADDRESS  REGISTERS 


ARO  -  The  increment  code  is  initialized  to  1.  The  address  portion  is 
initialized  to  the  address  of  the  AMP  buffer. 

AR1  -  The  increment  code  is  initialized  to  1.  The  address  portion  is 
initialized  to  the  address  of  the  VAR  buffer. 

AR2  -  The  increment  code  is  initialized  to  1.  The  address  portion  is 
initialized  to  the  address  of  the  RN  buffer. 

AR3  -  The  increment  code  is  initialized  to  1.  The  address  portion  is 
initialized  to  the  address  of  the  output  BB  buffer.  An  initial 
dummy  read  with  a  "subtract  increment"  directive  effectively 
initializes  the  address  portion  to  BB-1. 


f 
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INCREMENT  REGISTERS 

INC  0,  1,  2,  3  -  Wrap  code  must  be  1024. 

BRANCH  REGISTERS 

BR4  -  Unconditional  branch  to  return  to  SUPER  MACRO  at  completion  of 
BBAMP.  If  BBAMP  and  DSCLN  are  not  combined  by  a  SUPER  MACRO, 

BR4  should  branch  to  STOP. 

BCRO  -  Conditional  branch  over  the  variator  change  portion  of  the  main 
loop.  Conditioned  on  the  sign  of  the  incremented  count  CT. 

BNZ8  -  Branch  register  controlling  the  number  of  times  the  inner  loop  is 
executed. 

BNZ9  -  Branch  register  controlling  the  number  of  times  the  outer  loop  is 
executed.  Total  number  of  data  pairs  processed  is  product  of  inner 
times  outer  loop  executions.  BNZ9  allows  processing  of  arrays 
longer  than  512  points. 

SINE/COSINE  DESTINATION  REGISTERS 

Not  applicable. 

DECIMATE  REGISTER 

Not  applicable. 

SCALE  FACTOR  REGISTERS 

SFOH  -  All  True/Direct. 

SFOL  -  Left,  center,  and  right  prescalers  -  True/Direct,  postscaler  - 
True/Left  4. 

ARITHMETIC  ELEMENT  REGISTER  MAP 
See  Table  2-3. 

WORKING  STORE  MAP 
See  Table  2-4. 

PROGRAM  CODING  CHART 
See  Table  2-5. 
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BBAMP  AP  MICROPROGRAM  SCALING 


ADDER  IN 


MULTIPLIER  IN 


OUTPUT 
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TABLE  2-2  BBAMP  AEC  REGISTER  MAP 


Incr  Code  Address  Register  WS  Wrap  Increment  Register  Memory  Address  Register 

034  15  0  3  4  15 


1 


A(AMP) 


A(VF) 


0 

0 

1 

A(RN) 

1 

A(BB) 

BR/BCR  Unconditional/ 
Conditional  Branch 

Branch  Address 
0  15  ( 


BBAMP 


Sin/Cos  Destination  Register 
AEO  AE1  AE  2  AE  3 


78 


1516  21  22 


31 


Decimate  Register 
Count  Reset 


Comments 

N  =  number  of  spectral 
cells  to  process. 


BNZ  Branch  and  Count  Registers 
Count  Reset  Spare  Branch  Address 


2-4  BBAMP/DSCLN  WORKING  STORE  MAP 
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CHAPTER  3 

DISCRETE  LINE  SPECTRA  GENERATION  (DSCLN)  MACRO 

INTRODUCTION 

The  discrete  line  spectra  generation  algorithm  is  described  in  detail  in  a 
previous  publication.3  Briefly,  for  each  discrete  line  to  be  added  to  the  spectrum 
of  each  2048  sample  time  waveform  generated,  complex  components  are  added  to  six 
frequency  cells  -  three  on  each  side  of  the  desired  frequency.  The  magnitudes  and 
phases  of  the  six  complex  components  are  computed  from  the  desired  frequency  and 
amplitude  for  the  segment  to  be  generated  and  from  the  frequency  and  phase  of  the 
preceding  segment. 

The  exact  amplitude  weighting  for  the  six  components  is  given  by  the  periodic 
Hanning  function.4  For  DSCLN  MACRO,  a  three-segment  approximation  is  used  based 
on  sine  and  cosine  functions  derived  from  the  fractional  part  of  the  desired 
frequency.  Figure  3-1  illustrates  the  use  of  quarter-cycle  and  half-cycle  trig 
functions  to  generate  the  desired  weighting  function.  In  Figure  3-2,  for  a  desired 
spectral  line  of  amplitude  A  at  C  +  D  (C,  the  next  lower  cell  index  and  D,  the 
fractional  distance  to  the  specified  line  position)  the  following  computations  are 
performed  to  obtain  the  magnitudes  for  the  six  cell  components. 

For  cell  c  -  2: 

Mc_2  =  A( .035*sin(irD)*  cos  (4^)) 

=  A(.07*sin(^)*  cos2  (^) ) 

For  cell  c  -  1 : 

M(;_1  =  A(  .5  -  .5*sin(^)  +  .026*sin(uD)) 

=  A( .  5  -  .5*sin(4p)  +  .052*sin(^-)*cos(^)) 

3DAVIS,  R.  H.,  "Synthesis  of  Steady-State  Signal  Components  by  an  All-Digital 
System",  NOLTR  74-215,  Naval  Ordnance  Laboratory  (now  Naval  Surface  Weapons 
Center),  05  Dec  1974. 

4  Ibid. 
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CONSTRUCTION  OF  HANNING  APPROXIMATION 


FIGURE  3-2  DSCLN  ILLUSTRATION  OF  TERMINOLOGY 
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For  cell  c: 

Mc  =  A( .5  +  .  5*cos(^p)) 

For  cel  1  c  +  1 : 

Mc+]  =  A( . 5  +  .5*sin(^)) 

For  cell  c  +  2: 

Mc+2  =  A(’5  '  ■5co^f)  +  -O52*sin(^)*cos(^)) 

And  for  cell  c  +  3: 

Mc+3  =  A(.035*sinUD)*cos(^)) 

=  A(.07(sin2(^)*cos(^)). 

Given  D,  the  sin(-y-)  and  cos  are  computed  and  combined  to  form  the  above 
functions.  For  a  desired  phase  P,  the  component  phases  are  set  to  -P,  -P,  P, 

-P,  P,  P. 

Figure  3-3  is  a  BASIC  generated  plot  comparing  the  three-segment  approximation 
to  the  exact  Hanning  Function. 

FUNCTIONAL  DESCRIPTION 

DSCLN  MACRO  consists  of  two  phases  of  computation  for  each  discrete  line. 

The  first  phase  computes  the  center  frequency,  adds  the  FM  perturbation  to  obtain 
the  instantaneous  frequency,  and  computes  the  initial  phase  to  be  used  for  the  next 
iteration  of  DSCLN.  The  first  phase  also  computes  the  modulated  amplitude  to  be 
used  during  the  second  phase  computations  of  the  current  iteration.  The  second 
phase  computes  the  six  complex  spectral  components  for  each  line  based  on  the 
instantaneous  frequency  and  phase  stored  in  the  Discrete  Line  Table  (DLT)  input 
data  on  the  amplitude  just  obtained  during  Phase  1.  The  six  components  are  then 
added  to  the  appropriate  cells  of  the  spectrum  being  constructed. 


MATHEMATICAL  DESCRIPTION 

For  the  i^*1  line  during  the  j^  iteration,  the  following  computations  are 
performed: 

PhaseJj  CDOi  ^  j+2  =  CDOi  ^  j+1  +  SLi 

A.  L. 

where  CDO..  .  is  the  center  frequency  for  the  jtn  iteration  and  SL.  is  the  slew 

1  *  J  th  1 

rate  for  the  i  line.  Note  that  this  computation  produces  the  center  frequency 

CD0.j  valid  for  the  second  following  iteration 


FSM. 


j+1 


F$.  *  FM. 

»  '  9 


j+1 
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3-3  COMPARISON  OF  APPROXIMATION  AND  EXACT  HANNING  FUNCTION 
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where  FSMi  >  ^  is  the  appropriately  scaled  FM  perturbation  for  the  i 
the  (j+l)st  iteration. 


•  th 


line  during 


where  CD.  -  is  the  instantaneous  frequency  for  the  iteration, 
i  >  J 


where  A.  is  the  unmodulated  line  amplitude  for  the  i  line,  and  AM^ 


1 


•  th 


,th 


J 


is  the 


amplitude  modulation  factor  for  the  i  line  during  the  j  iteration. 


C°DPi,  -  [CD0j>  j+,  «■  FSM,_  .t,] 


Integer  part  plus  fraction  part  divided  by  4  -  where  the  integer  part  (H.O.  16  bits) 
is  the  center  cell  number  C  and  the  fractional  part  (L.O.  16  bits)  is  used  as  an 
angle  argument  <j>D  (in  BAM)  in  computing  the  Hanning  coefficients. 

Phase  2: 


Let  A  =  Ai  j 

K  =  cos(2tt  *  P) 


Z  =  sin(2ir  *  P) 

C  =  COs(2-rr  *  4>D) 

S  =  sin(2Tt  *  <j>D) 

7  =  .07 
5  =  .052 

BM2  (R,I)  =  Initial  conten 

=  Initial  conten 

=  Initial  conten 

=  Initial  conten 

=  Initial  conten 

=  Initial  conten 

=  Modified  contents  of  cell  C-2 
=  Modified  contents  of  cell  C-l 
=  Modified  contents  of  cell  C 
=  Modified  contents  of  cell  C+l 
=  Modified  contents  of  cell  C+2 
=  Modified  contents  of  cell  C+3 
The  following  computations  are  performed: 

M2R  =  BMf“  . . 

M2 1  =  bm; 


BM1  (R,I) 
BO  (R, I ) 
BP1  (R, I ) 
BP2  (R,I) 
BP3  (R,I) 
M2  (R,I) 
Ml  (R,I) 

B  (R,I) 

PI  (R.I) 
P2  (R.I) 
P3  (R.I) 


of 

cell 

O 

1 

ro 

of 

cell 

C-l 

of 

cell 

c 

of 

cell 

C+l 

of 

cell 

C+2 

of 

cell 

C+3 

MIR  =  BM1R  +  1 /2 { ASK-AK)  -  5SCAK 

Mil  =  BM1I  +  1 / 2 ( ASZ-AZ )  -  5SCAZ 

OR  =  BOR  +  1/2(AK+ACK) 

01  =  BOI  +  1 / 2 ( AZ+ACZ ) 

PIR  =  BP1R  -  1/2(AK+ASK) 

PI  I  =  BP1I  -  1/2 ( AZ+ASZ ) 

P2R  =  BP2R  +  1/2(ACK-AK)  +  5SCAK 
P2I  =  BP2I  +  1/2(ACZ-AZ)  +  5SCAZ 
P3R  =  BP3R  +  7CSSAK 

P3I  =  BP3I  +  7CSSAZ 

IMPLEMENTATION 

DSCLN  is  implemented  with  a  48-instruction  main  loop  which  processes  one  set 
of  Minor  Frame  (MNF)  computations  for  one  line.  An  eighteen-instruction  preamble 
generates  the  parameters  7  and  5  and  two  auxilaries,  Ml  and  PI,  used  in  the  main 
loop,  and  initializes  the  pipe  for  the  computations  for  the  first  line.  The  loop 
count  is  controlled  by  BNZ.A,  which  should  be  set  to  the  number  of  lines  to  be 
processed  minus  one.  Since  DSCLN  is  the  last  AP  program  to  be  executed  under  the 
SUPER  MACRO,  SMI,  a  halt  is  executed  after  completion. 

Since  intermediate  results  CDO,  CODP  and  P  are  updated  in  Working  Store  by 
DSCLN,  this  area  of  WS  must  be  returned  to  Bulk  Store  upon  completion  of  DSCLN. 

SCALING 

See  Table  3-1 . 


ARITHMETIC  ELEMENT  CONTROLLER  (AEC)  REGISTERS 
See  Table  3-2. 

ADDRESS  REGISTERS 

ARO  -  The  initial  value  is  set  to  cos"'' (7)  =  27C8.  During  main  loop 

computations  ARO  contains  the  spectral  cell  addresses  C-2,  ...,  C+3. 

ART  -  The  initial  value  is  set  to  sin"' (5)  =  2809.  During  main  loop 

computations  AR1  contains  the  fraction  $D  from  which  the  coefficients 
C  and  S  are  obtained. 

AR2  -  During  main  loop  computations  AR2  is  loaded  with  FMA  and  used  to 
access  the  MTR  functions  for  FM  values. 
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AR3  -  Ouring  main  loop  computations  AR3  is  loaded  with  AMA  and  used  to 
access  the  MTR  functions  for  AM  values. 

AR4  -  The  increment  code  is  initialized  to  1.  The  address  portion  is 
initialized  to  the  starting  address  of  the  SLP,  CDO  data  area. 

AR5  -  The  increment  code  is  initialized  to  1.  The  address  portion  is 
initialized  to  the  starting  address  of  the  FAMA,  FSA  data  area. 

AR6  -  The  increment  code  is  initialized  to  1.  The  address  portion  is 
initialized  to  the  starting  address  of  the  CODP  data  aiea. 

INCREMENT  REGISTERS 

INCRO  -  1024  Wrap.  The  increment  value  is  set  to  -5. 

INCR1  -  1024  Wrap.  The  increment  value  is  set  to  -2. 

INCR2  -  1024  Wrap.  The  increment  value  is  set  to  +2. 

INCR3  -  1024  Wrap. 

INCR4  -  1024  Wrap.  The  increment  value  is  set  to  +2. 

INCR5  -  1024  Wrap. 

INCR6  -  1024  Wrap. 

BRANCH  REGISTERS 

BNZA  -  Branch  and  Count  Register  controlling  the  number  of  executions  of 
the  main  loop.  Count  and  Reset  fields  are  initialized  to  N-l , 
where  N  is  the  number  of  lines  to  be  processed.  The  Branch  Address 
is  DSCLP. 

SINE/COSINE  DESTINATION  REGISTERS 

SCDO  -  Destine  Cosine  only  to  AEO. 

SCDl  -  Destine  Cosine/Sine  to  AEO. 

SCD2  -  No-op  (00). 

SCD3  -  Destine  Cosine/Sine  to  AEO. 

DECIMATE  REGISTER 

Not  Applicable. 

SCALE  FACTOR  REGISTERS 

SFOH  -  All  True-Direct. 

SFOL  -  Left  and  center  prescalers  -  True/Right  1,  right  prescaler  and 
postscaler  -  True/Direct. 
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SF1H  -  Left  prescaler  -  True/Right  4,  center  and  right  prescalers  - 
True/Right  2,  postscaler  True/Left  2. 

SF1L  -  All  prescalers  -  True/Direct,  postscaler  -  True/Left  4. 

SF2H  -  Left  prescaler  -  True/Right  1,  center  and  right  prescalers  - 
True/Direct,  postscaler  -  True/Left  1. 

SF2L  -  All  prescalers  -  True/Direct,  postscaler  -  True/Left  2. 

SF3H  -  Left  prescaler  -  True/Right  2,  center  prescaler  -  True/Direct, 
right  prescaler  -  True/Right  1,  postscaler  -  True/Left  1. 

SF3L  -  Left  and  center  prescalers  -  True/Right  2,  right  prescaler  - 
True/Direct,  postscaler  -  True/Left  2. 

ARITHMETIC  ELEMENT  REGISTER  MAP 

See  Table  3-3. 

WORKING  STORE  MAP 

See  Table  2-4.  (BBAMP  and  DSCLN  are  combined.) 

PROGRAM  CODING  CHART 
See  Table  3-4. 
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TABLE  3-1  OSCLN  AP  MICROPROGRAM  SCALING 


LABEL 

ADDER  IN 

MULTIPLIER  IN 

OUTPUT 

INPUTS: 

CO/ DP 

HO,  HI 4 

CDO 

FO 

SL 

H16  (L.O.HW) 

P 

H16  (L.O.HW) 

A 

HI  2 

FM 

HI  2 

AM 

HI  5 

" 

FS 

H8 

7 

HI  7 

5 

HI  8 

8M2R,  BM2I 

H12 

! 

BM1R,  SMI  I 

HI  2 

BOR,  BOI 

HI  2 

BP1R,  BP1 I 

HI  2 

BP2R,  BP2I 

H12 

BP3R,  BP3I 

HI  2 

INTERMEDIATE  RESULTS: 

S,  C 

H14 

Z,  K 

H 1 4 

7,  5 

H18 

A 

HI  2 

7C 

HI  6 

5S 

Hi  6 

CS 

HI  4 

3-10 


AZC _ 


7CAK 


7CAZ 


5SCAK 


5SCAZ 


7CCSAK 


7CCSAZ 


7SSCAK 


7SSCAZ 


Tl,  T2 


T3,  T4 


FSM 


D4 


SL 


PP 


NCD4/NCD 


CDP4 


NC 


cc 
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TABLE  3-1 

(CONTINUED! 

LABEL 

ADDER  IN 

MULTIPLIER  IN 

OUTPUT 

PSMd 

1  FO 

1 

1 

OUTPUTS : 

M2R.  M2 1 

H12 

MIR.  Mil 

HI  2 

OR.  01 

H12 

P1R.  P1I 

HI  2 

P2R,  P2I 

HI  2 

P3R,  P3I 

H12 
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TABLE  3-2  DSCLN  AEC  REGISTER  MAP 

Incr  Code  Address  Register  WS  Wrap  Increment  Register  Memory  Address  Register 

3 


27C8=Cos  (7) 


2809=Sin'1(5) 


0  FFB=-5 


0  FFE=-2 


Sin/Cos  Destination  Register 
AEO  AE1  AE2  AE  3 


1 

620=A(SL,P,CD0) 

1 

6E8=A(FAMA,FSA) 

1 

58C=A(C0OP) 

0 

002 

0 

0 

0 

0 

— 

BR/BCR  Unconditional/ 

Conditional  Branch  BN2  Branch  and  Count  Registers 


Branch  Address 
0 _  15 


Count  Reset  Spare  8ranch  Address 
0  78  1516  21  22  31 


1031 


Decimate  Register 
Count  Reset 


Comments 

N  =  *  of  lines 


1\ 
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TABLE  3-4  (CONTINUED) 
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CHAPTER  4 

SINGLE  REAL  FFT  (SRFFT)  MACRO 


FUNCTIONAL  DESCRIPTION 


The  SRFFT  MACRO  performs  the  complex  multiplications  required  to  allow  the 
FFT  to  effect  the  transformation  from  real  time  series  to  single-sided  spectra. 
SRFFT  is  implemented  as  a  separate  pass  which  must  be  executed  on  a  single-sided 
spectrum  before  bit  reversal ,  or  as  a  last  step  after  a  forward  transform  and  bit 
reversal  to  produce  a  single-sided  spectrum.  For  N  =  1024,  2048  time  samples 
are  produced  from  1024  spectral  samples  or  vice-versa,  and  the  six  passes  (1  SRFFT, 
5  FFT4)  require  2.46  ms. 

The  computation  involves  pairs  of  complex  samples  from  symmetrical  locations 

in  the  data  array,  i.  e.,  outputs  Fp  +  FN_n  are  derived  from  inputs  Gp  and  GN_n- 

l'2-trr 
2N 


The  computation  requires  the  exponential  e  ,  which  is  generated  by  the  Sin/Cos 


Generator. 


MATHEMATICAL  DESCRIPTION 

The  input  array  is  denoted  by  Gn>  n=0,  ...,  N-l.  The  output  array  is  denoted 
by  Fn,  n=0,  N-l.  The  computations  are  broken  into  the  following  steps,  where 
GRi ,  GI.,  FRi ,  FI..,  etc.,  denote  the  real  and  imaginary  parts  of  Gi  +  F^ , 
respectively,  and  C.  and  S^  the  real  and  imaginary  parts  of  the  complex  exponential 
exp  (j2tt  2ft): 

=  (GRn_.  -  GR.)/2 

CIi  =  (GIN-i  -  GV/2 

SCR.=  S.  *  CR. 

CCR.=  C.  *  CRi 

CCI.=  ci  *  Ch 
SCI^  s.  *  CI. 

ARn-  =  (GRn_.  +  GR.)/2 

AI i  =  (GIN-1  -  GI-5/2 
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FRi  =  (SCRi  +  CCK  +  ARi/4)*4 

FI,.  =  (SC I.  -  CCRi  +  AIi/4)*4 

FRN.i  =  '  CCIi  +  ARi/4)*4 

FIN_i  =  (SCI i  -  CCRi  -  AI./4)*4 
IMPLEMENTATION 

The  SRFFT  MACRO  consists  of  a  sixteen-instruction  prologue  and  an  eight- 
instruction  main  loop.  During  the  prologue,  two  minus-one's  are  created  and  placed 
in  MLR.2H  and  MRR.1H  to  allow  data  to  be  recycled  through  the  multiplier  with  only 
a  sign  change.  The  (N-i)  output  address  counter  AR.5  is  artificially  incremented  by 
one  to  allow  the  first  output  to  be  effected  in  the  main  loop  with  a  decrement 
associated  with  it.  Also  unique  to  the  prologue  is  the  computation  for  Fg,  which 
for  a  forward  transform  is  FRO  =  GRO/2  +  GIO/2,  and  for  an  inverse  transform  is 
FRq  =  GRq,  Fig  =  GRg.  These  are  selectable  by  SF2H  and  SF2L. 

Scaling  is  applied  in  the  computation  of  AR,  Al,  CR,  and  Cl  (SFOL)  to  effect 
the  overall  factor  of  1/2  that  is  required  from  input  to  output.  In  addition, 

SF1H  is  used  to  compensate  for  the  effective  scaling  index  difference  of  two  that  is 
introduced  when  CR  and  Cl  are  multiplied  by  C  and  S.  If  an  additional  net  scaling 
is  desired,  it  may  be  effected  by  reducing  the  number  of  postscaler  shifts  in 
SF1H. 

The  main  loop  count  control  is  based  on  BNZ8,  with  an  outside  loop  counter 
BNZ9  to  permit  operation  on  arrays  larger  than  512  points. 

SCALING 

See  Table  4-1 . 

ARITHMETIC  ELEMENT  CONTROLLER  (AEC)  REGISTERS 

See  Table  4-2. 

ADDRESS  REGISTERS 

ARO  -  The  initial  value  is  set  to  zero,  including  increment  code. 

AR2  -  The  increment  code  is  set  to  1.  The  address  portion  is  initialized 
to  the  starting  address  of  the  Input  Array. 

AR3  -  The  increment  code  is  set  to  1.  The  address  portion  is  initialized 
to  the  last  address  of  the  Input  Array. 

AR4  -  The  increment  code  is  set  to  1.  The  address  portion  is  initialized 
to  the  starting  address  of  the  Output  Array. 
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AR5  -  The  increment  code  is  set  to  1.  The  address  portion  is  initialized 
to  the  last  address  of  the  Output  Array. 

INCREMENT  REGISTERS 

INCO  -  Set  to  1/2N  as  BAM  angle  increment. 

INC2,  3,  4,  5  -  Wrap  code  must  be  set  to  1024. 

BRANCH  REGISTERS 

BNZB  -  Branch  and  Count  Register  controlling  the  number  of  executions  of 
the  inner  loop.  Count  and  Reset  fields  are  initialized  to  N/4-2 
and  N/4-1,  respectively,  where  N  is  the  array  size.  The  Branch 
Address  is  SRFLP. 

BNZC  -  Branch  and  Count  Register  controlling  the  number  of  executions  of 
the  outer  loop.  Count  and  Reset  fields  are  initialized  to  1.  The 
Branch  Address  is  SRFLP. 

SINE/COSINE  DESTINATION  REGISTERS 

SCDO  -  Destine  Sine/Cosine  to  AEO  and  AE1. 

DECIMATE  REGISTERS 

Not  applicable. 

SCALE  FACTOR  REGISTERS 

SFOH  -  All  True/Direct. 

SFOL  -  Left,  center  and  right  prescalers  -  True/Right  1,  postscaler  - 
True/Direct. 

SF1H  -  Left  prescaler  -  True/Direct,  center  prescaler  -  True/Right  2, 
right  prescaler  -  True/Direct,  postscaler  -  True/Left  2. 

SF2H  -  Forward  Transform:  Left  and  center  prescalers  -  True/Right  1, 
right  prescaler  and  postscaler  -  True/Direct. 

Inverse  Transform:  Left  and  center  prescalers  -  Inhibit,  right 
prescaler  and  postscaler  -  True/Direct. 

SF2L  -  Forward  Transform:  Left  prescaler  -  True/Direct,  center  and  right 
prescalers  -  Inhibit,  postscaler  -  True/Direct. 

Inverse  Transform:  Left  prescaler  -  True/Direct,  center  and  right 
prescalers  -  Inhibit,  postscaler  -  True/Direct. 
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ARITHMETIC  ELEMENT  REGISTER  MAP 
See  Table  4-3. 

WORKING  STORE  MAP 
See  Table  4-4. 


PROGRAM  CODING  CHART 
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TABLE  4-1  SRFFT  AP  MICROPROGRAM  SCALING 


LA8EL 

ADDER  IN 

MULTIPLIER  IN 

OUTPUT 

INPUTS: 

GR.  GI 

HI  2 

H12 

GRN,  GIN 

HI  2 

HI  2 

C,  S 

H1A 

INTERMEDIATE  RESULTS: 

MGRN,  MG  IN 

H 1 2  (L.O.HW) 

MGR,  MG I 

H12  (L.O.HW) 

CR,  Cl 

HI  2 

AR,  AI 

FI  2 

1 

CCR,  CCI 

1  FI  0 

SCR,  SCI 

FI  0 

OUTPUTS: 

■ 

FR,  FI 

! 

Hi  2 

FRN,  FIN 

HI  2 
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TABLE  4-2  SRFFT  AEC  REGISTER  MAP 


Incr  Code  Address  Register  US  Wrap  Increment  Register  Memory  Address  Registe 

034  15  0  3  4  15 


■ 

1 

nr 

■ 

■ 

A(  INPUT)  =  I 

■ 

A(EOI)  =  I+N-l 

1 

A (OUTPUT)  =  0 

1 

A(EOO)  =  O+N-l 

SRFFT 


H 

0 

Sin/Cos  Destination  Regis 
AEO  AE1  AE2  AE3 


no 

no 

N 

M 

m 

m 

H 

m 

E; 

HR 

m 

■ 

Decimate  Register 

0 

Count  Reset 

Comments 


8R/BCR  Unconditional/ 

Conditional  Branch  8NZ  Branch  and  Count  Registers 


N  -  Number  of  Complex 
Points 


Branch  Address 
0  15 


Count  Reset  Spare  Branch  Address 
0  78  1516  21  22 _ 31 


E01  =  End  of  Input 
EOO  *  End  of  Output 


-IwlMfe-aJArtL. 
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CHAPTER  5 

SCALE  AND  ADD  TWO  ARRAYS  (SCLA2)  MACRO 

FUNCTIONAL  DESCRIPTION 

SCLA2  scales  an  array  X  by  a  constant  A,  scales  an  array  Y  by  a  constant  B, 
and  combines  the  two  scaled  arrays  on  an  element-by-element  basis  to  form  an  output 
array  Z.  X,  Y,  and  Z  must  be  of  the  same  number  of  elements  N.  The  scaling 
coefficients  A  and  B  are  obtained  via  the  Sin/Cos  Generator  in  order  to  avoid  the 
requirement  for  two  words  of  Working  Store  (WS).  These  must  therefore  be  loaded 
from  the  corresponding  CFCB  as  the  ARCCOS  (A)  and  the  ARCCOS  (B),  respectively. 

MATHEMATICAL  DESCRIPTION 

Z.  =  A  *  Xi  +  B  *  Y.;  i  =  1,  2,  ....  N 

where  A  and  B  are  two  scalers,  and  Xi ,  Y,.  and  Z^  are  the  ith  elements  of  arrays 
X,  Y  and  Z,  respectively. 

IMPLEMENTATION 

SCLA2  is  implemented  with  a  six-instruction  main  loop  which  processes  one 
pair  of  X's  and  one  pair  of  Y's  to  form  one  pair  of  Z's.  A  twelve-instruction 
preamble  obtains  the  coefficients  A  and  B  via  the  Sin/Cos  Generator  and  loads 
the  pipe  with  the  first  pairs  of  X's  and  Y’s.  Each  pass  through  the  main  loop 
then  stores  one  pair  of  Z's  in  the  output  array.  Since  the  maximum  count  is  255, 
or  256  passes  through  the  loop,  an  outside  loop  is  required  to  permit  array  sizes 
N  greater  than  512  elements.  The  inside  loop  is  controlled  by  BNZ8  and  the 
outside  loop  is  controlled  by  BNZ9. 

SCALING 

See  Table  5-1 . 

ARITHMETIC  ELEMENT  CONTROLLER  (AEC)  REGISTERS 
See  Table  5-2. 
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ADDRESS  REGISTERS 

ARO  -  ARCCOS(A),  i.  e. ,  ARO  is  set  to  an  angle  a  such  that  cos(a)  =  A. 
AR1  -  ARCCQS(B) . 

AR2  -  The  increment  code  is  initialized  to  one.  The  address  is  in¬ 
itialized  to  the  first  address  of  the  X  array. 

AR3  -  The  increment  code  is  initialized  to  one.  The  address  is  in¬ 
itialized  to  the  first  address  of  the  Y  array. 

AR4  -  The  increment  code  is  initialized  to  one.  The  address  is  in¬ 
itialized  to  the  first  address  of  the  Z  array. 

INCREMENT  REGISTERS 

Not  applicable.  Wrap  codes  must  be  1024  for  INC  2,  3,  4. 

BRANCH  REGISTERS 

BNZ8  -  Branch  and  Count  register  controlling  the  inner  loop.  Reset  and 
Count  fields  are  initialized  to  N/M-l,  where  N  is  the  number  of 
elements  in  each  of  the  arrays,  and  M  is  the  number  of  times  the 
outside  loop  is  executed.  Branch  Address  is  SC2LP. 

BNZ9  -  Branch  and  Count  register  controlling  the  number  of  executions  of 
the  outer  loop.  Reset  and  Count  fields  are  initialized  to  M-l . 
Branch  Address  is  SC2LP. 

SINE/COSINE  DESTINATION  REGISTERS 

SCDO  -  Destine  Cosine  only  to  both  AEO  and  AE1. 

SCD1  -  Destine  Cosine  only  to  both  AEO  and  AE1 . 

DECIMATE  REGISTER 

Not  applicable. 

SCALE  FACTOR  REGISTERS 

SFOH  -  Left,  center  and  right  prescalers  -  True/Direct.  Postscaler  - 
True/ Left  2. 

ARITHMETIC  ELEMENT  REGISTER  MAP 
See  Table  5-3. 
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WORKING  STORE  MAP 
See  Table  5-4. 
PROGRAM  CODING  CHART 
See  Table  5-5. 
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CHAPTER  6 

ASYNCHRONOUS  SAMPLE,  SCALE  AND  SUM  (ASSSS)  MACRO 


FUNCTIONAL  DESCRIPTION 

The  ASSSS  MACRO  performs  the  dopplered  resampling  of  the  input  signal  D  for 
each  multipath  arrival  at  the  sensor.  Linear  interpolation  is  employed  to  provide 
some  protection  against  aliasing,  and  is  equivalent  to  passing  the  signal  through 
a  filter  with  a  frequency  response  of: 


H(f )  = 


.2/  Ttf  \ 

sin  (-F— ) 


(#■>' 


The  parameter  f$  is  the  sample  rate  of  the  input  or  source 


before  resampling. 

signal.  The  resampled  signal  is  multiplied  by  a  gain  parameter  A  and  added  to  a 
Sum  array  S.  Since  there  is  not  a  one-to-one  correspondence  between  the  samples 
of  the  input  signal  D  and  the  output,  to  add  N  new  samples  to  the  Sum  array  may 
require  more  or  less  than  N  samples  of  D.  The  resampling  may  start  at  an  arbitrary 
phase  (or  fractional  position  F,  0  <_  F  <  1)  with  respect  to  the  input  samples,  but 
the  MACRO  is  intended  for  application  only  where  the  resample  rate  is  relatively 
close  to  the  original  rate,  say  +  10%  maximum. 

MATHEMATICAL  DESCRIPTION 

Phase  1 :  Apply  amplitude  factor  to  source  input  samples 
AH  =  A  *  OH 
AL  =  A  *  DL 

where  A  is  the  amplitude  factor  and  DH  and  DL  are  the  input  samples  in  the  high 
and  low  halves  of  the  Working  Store  (WS)  word. 

Phase  2:  Form  interpolation  fractions  for  next  pair  of  outputs. 

FI  =  F2  +  D  +  0C 
F3  =  -F2  -  D  +  1C 

F2  =  FI  +  0  +  0C 
F4  =  -FI  -  D  +  1C 
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where  D  is  the  resample  rate  parameter,  FI  and  F2  are  the  resample  fractions,  F3 
and  F4  are  the  complements  of  FI  and  F2,  respectively,  with  respect  to  unity,  and 
OC  and  1C  are  conditional  "ones"  that  are  used  to  constrain  the  F's  to  interval  0-1. 

Phase  3:  Form  the  products  of  the  scaled  inputs  with  the  interpolation 
fractions. 

1AH  =  FI  *  AH 
2AL  =  F2  *  AL 
3AL  =  F3  *  AL 
4 AH  =  F4  *  AH 

Phase  4:  Combine  outputs  of  Phase  3  with  sample  values  from  Sum  array  to  form 
final  results. 

For  Inphase  loop: 

OH  =  SH  +  3AL  +  1  AH 
0L  =  SL  +  4AH  +  2AL 
For  Outphase  loop: 

OH  =  SH  +  4AH  +  2AL 
0L  =  SL  +  3AL  +  1AH 

IMPLEMENTATION 

The  ASSSS  MACRO  consists  of  two  main  loops,  each  eight  instructions  long, 

eight  transition  sequences  for  the  eight  possible  ways  that  transitions  between 

the  two  major  loops  can  occur,  plus  some  initialization  and  finalization  code. 

If  the  block  size  is  1024  samples  (512  words),  both  the  input  array  and  the  output 

array  may  exceed  512  words,  and  the  MACRO  must  operate  from  one  half  of  Working 

Store  to  the  other  half.  Smaller  block  sizes  may  be  handled  within  one  side  of  WS. 

The  two  main  loops  are  designated  the  "Inphase"  loop  and  the  "Outphase"  loop, 

depending  on  how  the  input  array  and  output  array  are  related.  Actually,  a  more 

reasonable  nomenclature  mighx  be  "leading"  and  "lagging"  loops,  since  for  the 

"inphase"  or  "leading"  loop,  the  two  samples  in  an  output  data  word  are  derived 

from  the  two  samples  in  the  input  data  word  plus  one  contribution  from  the  previous 

input  data  word.  For  the  "outphase"  or  "lagging"  loop,  the  output  data  samples 
» 

are  derived  from  the  two  samples  in  the  input  data  word  plus  one  contribution  from 
the  following  input  data  word. 


I 
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Since  the  transitions  will  occur  at  rates  less  than  one  in  twenty  or  more 
times  through  the  main  loops,  the  timing  is  nominally  determined  by  the  eight  steps 
per  word  for  each  of  the  main  loops.  For  nominally  512  executions  of  one  or  the 
other  of  the  main  loops  (1024  samples),  the  approximate  execution  time  for  the 
ASSSS  MACRO  is  .41  ms. 

The  operation  of  the  two  main  loops  is  illustrated  in  Figures  6-1  and  6-2. 

Each  figure  follows  the  two  samples  in  one  input  word  through  the  pipeline  to  the 
three  output  samples  that  they  affect.  Since  the  outputs  contain  terms  that  involve 
triple  products  of  an  interpolation  fraction,  an  amplitude  and  the  input  data  sample, 
two  full  trips  through  the  pipeline  are  required,  beginning  with  a  read  of  the 
input  word  in  the  first  trip  through  the  loop  and  ending  with  the  write  of  the  last 
affected  output  in  the  fourth  following  pass  through  the  loop.  On  the  first  pass 
through  the  pipeline,  the  input  data  samples  are  multiplied  by  the  amplitude  factor. 
Approximately  concurrent  with  these  first  products  being  passed  through  the  adder, 
the  four  appropriate  interpolation  fractions  are  generated  and  passed  to  the 
multiplier.  On  the  second  pass  through  the  pipeline,  the  scaled  samples  are 
multiplied  by  the  interpolation  fractions  and  these  triple  products  are  combined 
with  the  corresponding  Sum  array  samples  to  form  the  final  composite  output. 

At  the  time  the  interpolation  fractions  F.  are  formed,  a  test  is  made  to 
determine  whether  they  have  progressed  out  of  the  range  0  <  F.  <  1.  Depending  on 
which  of  the  two  fractions  FI  or  F2  have  exceeded  the  valid  range,  and  whether  the 
increment  D  is  positive  or  negative,  the  program  deviates  through  one  of  eight 
transition  routines  to  resume  operation  in  the  alternate  (inphase  or  outphase)  loop. 
These  transition  routines  re-adjust  the  appropriate  interpolation  fractions  to  the 
valid  range  and  compute  the  one  or  two  output  values  that  do  not  conform  to  the 
procedures  in  the  main  loops.  The  following  gives  a  brief  description  of  each  of 
the  transition  sequences. 

PIQ01  In  to  Out  On  FI ,  D  positive 

At  the  time  the  branch  is  taken,  the  following  quantities  have  been  erroneously 
computed: 

FI  =  -1  +  e  F3  =  1  -  e  -  2D 

F2=e+D  F4  =  1  -  e  -  D 

where  0  is  the  amount  the  fraction  F  is  incremented  each  sample,  and  e  is  a  small 
number  between  0  and  D.  The  following  two  transition  outputs  are  computed: 

OH  *  SH  =  e  *  AL  +  (1-e)  *  AH  (AH  from  current  input) 
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and  OL  =  SL  +  (e+D)  *  AH  +  (1-e-D)  *  AL  (AH  from  next  input) 

which,  using  the  fractional  quantities  already  computed,  become 
OH  =  SH  +  FI  *  AL  -  FI  *  AH  +  AL 

and  OL  =  SL  -  F2  *  AH  -  F4  *  AL 

Finally,  the  next  interpolation  fractions  are  computed 
F2  =  e  +  2D  F4  =  1  -  e  -  2D 

and  computation  resumes  in  the  Outphase  loop  at  M0PL2. 

PI 002  In  to  Out  On  F2 

At  the  time  the  branch  is  taken,  the  following  quantities  have  been  erroneously 
computed: 

F2  =  -1  +  e  F4  =  -e 

FI  =  e  +  D  F3  =  1  -  e  -  D 

The  transitional  outputs  become 

OL  =  SL  +  e  *  AH  +  (1-e)  *  AL  (previous  AL) 

or  OL  =  SL  -  F4  *  AH  -  F2  *  AL 

and  OH  =  SH  +  (e+D)  *  AL  +  (1-e-D)  *  AH  (current  AL) 

or  OH  =  SH  +  FI  *  AL  +  F3  *  AH 

The  adjusted  interpolation  fractions  are: 

F2  =  e  +  2D  F4  =  1  -  e  -  2D 

Computation  resumes  at  M0PL1  in  the  Outphase  loop. 

P0I01  Out  to  In  On  FI 

At  the  time  the  branch  is  taken,  the  following  quantities  have  been  erroneously 
computed: 

FI  =  -1  +  e  F3  =  -e 

F2=e+D  F4  =  1  -  e  -  D 

The  transition  output  is 

OL  =  SL  +  e  *  AL  +  (1-e)  *  AH 
or  OL  =  SL  -  F3  *  AL  -  FI  *  AH 
and  the  adjusted  interpolation  fractions  are 

FI  =  e  +  D  F3  =  1  -  e  -  D 

Computation  resumes  at  MIPL1  in  the  Inphase  loop. 

PQI02  Out  to  In  On  F2 
Erroneously  computed  fractions  are: 

F2  =  -1  +  e  F4  =  1  -  e  -  2D 

FI  =  e  +  D  F3  =  1  -  e  -  D 
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(AL  from  previous  input) 


2D 


The  transition  output  is: 

OH  =  SH  -  e  *  AH  +  (1-e)  *  AL 
or  OH  =  SH  +  F2  *  AH  +  AH  -  F2  *  AL 
The  adjusted  interpolation  fractions  are: 

FI  =  e  +  2D  F3  =  1  -  e 

Computations  resume  at  MIPL2  in  the  Inphase  loop. 

The  four  transition  sequences  for  negative  fraction  increment  D  all  produce 
transition  outputs  which  are  not  exact.  Each  of  the  outputs  is  in  error  by  a  term 
involving  a  fraction  factor  that  has  a  magnitude  of  (-2D)  or  less.  Since  by 
definition  D  is  much  less  than  unity,  this  term  can  safely  be  ignored.  The 
alternative  would  have  required  long  transition  sequences  involving  "backing  up" 
the  computations  to  obtain  data  which  is  no  longer  in  the  pipeline. 

MI0Q1  In  to  Out  On  FI ,  D  negative 

At  the  time  the  branch  is  taken,  the  following  quantities  have  been  erroneously 
computed: 

FI  =  -e  .-2  =  1  -  e  +  D 

F3  =  -1  +  e  F4  -  e  -  D 

The  transition  outputs  are: 

OH  =  SH  -  F3  *  AL 
=  SH  +  (1-e)  *  AL 

(OH  =  SH  +  (1-e)  *  AL  +  e  *  AH  is  exact) 
and 

OL  =  SL  +  F2  *  AH 
=  SL  +  (1-e+D)  *  AH 

(OL  =  SL  +  (1-e+D)  *  AH  +  (e-D)  *  AL  is  exact). 

The  adjusted  interpolation  fractions  are: 

F4  =  e  -  2D  F3  =  e  -  3D 

F2  =  1  -  e  +  2D 

Computation  resumes  at  M0PL1  in  the  Outphase  loop. 

MI002  In  to  Out  On  F2,  D  negative 

At  the  time  the  branch  is  taken,  the  following  quantities  have  been  erroneously 
computed: 

T2  =  -e  F4  =  -1  +  e 

FI  =  1  -  e  +  D  F3  =  e  -  D 

The  transition  output  is: 
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OL  =  SL  -  F4  *  AH 
=  SL  +  (1-e)  *  AH 

(OL  =  SL  +  (1-e)  *  AH  +  e  *  AL  is  exact). 

The  adjusted  interpolation  fractions  are: 

F2  =  1  -  e  +  D  F3  =  e  -  2D 

Computation  resumes  at  M0PL1  in  the  Outphase  loop. 

MOIOl  Out  to  In  On  FI ,  D  negative 

At  the  time  the  branch  is  taken  the  following  quantities  have  been  erroneously 
computed: 

FI  =  -e  F3  =  -1  +  e 

F2  =  1  -  e  +  D  F4  =  e  -  D 

The  transition  output  is: 

OL  =  SL  -  F3  *  AL 
=  SL  +  (1-e)  *  AL 

(OL  =  SL  +  (1-e)  *  AL  +  e  *  AH  is  exact). 

The  adjusted  interpolation  fractions  are: 

F4  =  e  -  2D 
FI  =  1  -  e  +  0 

Computation  resumes  at  MIPL1  in  the  Inphase  loop. 

MI 002  Out  to  In  On  F2,  D  negative 

At  the  time  the  branch  is  taken  the  following  quantities  have  been  erroneously 
computed: 

F4  =  -1  +  e  F2  =  -e 

F3  =  e  -  D  FI  »  1  -  e  +  0 

F4  =  e  -  2D 

The  transition  outputs  are: 

OH  =  SH  -  4AH 

=  SH  +  (1-e)  *  AH 

(OH  =  SH  +  (1-e)  *  AH  +  e  *  AL  is  exact) 
and  OL  =  SL  +  1 AL 

=  SL  +  (1-e+D)  *  AL 

(OL  =  SL  +  (1-e+D)  *  AL  +  (e-D)  *  AH  is  exact). 

The  adjusted  interpolation  fractions  are: 

FI  =  1  -  e  +  2D 
F4  =  e  -  3D 
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Computation  resumes  at  MIPLl  in  the  Inphase  loop. 
SCALING 

See  Table  6-1 . 


ARITHMETIC  ELEMENT  CONTROLLER  (AEC)  REGISTERS 


See  Table  6-2. 


ADORESS  REGISTERS 

ARO  -  The  increment  code  is  set  to  one.  The  address  field  is  initialized 
to  the  location  of  the  gain  parameter  A. 

ART  -  The  increment  code  is  set  to  one.  The  address  field  is  initialized 
to  the  location  of  the  F  and  D  parameters. 

AR2  -  The  increment  code  is  set  to  one.  The  address  field  is  initialized 
to  the  location  of  the  resampled  input. 

AR3  -  The  increment  code  is  set  to  one.  The  address  field  is  initialized 
to  the  location  of  the  sum  input. 

AR4  -  The  increment  code  is  set  to  one.  The  address  field  is  initialized 
to  the  location  of  the  resampled  output. 

INCREMENT  REGISTERS 

Not  applicable.  Wrap  codes  must  be  1024  for  INC  2,  3,  4. 

BRANCH  REGISTERS 

BCRO  -  Conditional  branch.  For  D  positive,  branch  to  PI001  on  FI.  For 
D  negative,  branch  to  MI001  on  FI. 

BCR1  -  Conditional  branch.  For  D  positive,  branch  to  PI002  on  F2.  For 
D  negative,  branch  to  MI002  on  F2. 

BR2  -  Unconditional  branch  to  MIPL2. 

BR3  -  Unconditional  branch  to  M0PL2. 

BCR4  -  Conditional  branch.  For  D  positive,  branch  to  P0I01  on  FI.  For 
D  negative,  branch  to  MO 101  on  FI. 

BCR5  -  Conditional  branch.  For  D  positive,  branch  to  P0I02  on  F2.  For 
0  negative,  branch  to  M0I02  on  F2. 

BR6  -  Unconditional  branch  to  MIPLl. 

BR7  -  Unconditional  branch  to  M0PL1 . 

BNZ8  -  Branch  and  Count  register  controlling  the  number  of  executions  of 
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the  Inphase  processing  loop.  The  Branch  Address  is  MIPLP. 

BNZ9  -  Branch  and  Count  register  controlling  the  number  of  executions  of 
the  Outphase  processing  loop.  The  Branch  Address  is  MOPLP. 

BNZA  -  Unconditional  branch  to  MIPLO. 

BNZB  -  Unconditional  branch  to  MOPLO. 

BNZC  -  Unconditional  branch  to  STOP. 

SINE/COSINE  DESTINATION  REGISTERS 

Not  applicable. 

DECIMATE  REGISTER 

Not  applicable. 

SCALE  FACTOR  REGISTERS 

SFOH  -  All  True/Direct. 

SFOL  -  All  prescalers  -  True/Direct.  Postscaler  -  True/Left  1. 

SF1H  -  Left  prescaler  -  True/Right  4;  center  and  right  prescalers  - 
True/Direct;  postscaler  -  True/Left  4. 

ARITHMETIC  ELEMENT  REGISTER  MAP 

See  Table  6-3. 

WORKING  STORE  MAP 

See  Table  6-4. 

PROGRAM  CODING  CHART 

See  Table  6-5. 
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TABLE  6-1  ASSSS  AP  MICROPROGRAM  SCALING 


LABEL 

ADDER  IN 

MULTIPLIER  IN 

OUTPUT 

INPUTS : 

A 

HI  2 

F 

FI  5 

D 

FI  5 

DH,  DL 

HI  2 

SH,  SL 

H12 

INTERMEDIATE  RESULTS: 

AH,  AL 

F8 

HI  2 

FI,  F2 

FI  5 

HI  5 

F3,  F4 

FI  5 

Hi  5 

1AH,  2AL 

m 

3 AH,  4AL 

FI  1 

OT 

FI  2 

1 AL,  2AH 

FI  1 

OUTPUTS: 

OH,  OL 

HI  2 
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TA8LE  6-2  ASSSS  AEC  REGISTER  MAP 
Incr  Code  Address  Register  MS  Wrap  Increment  Register 


Memory  Address  Register 


0  3  4 


15  0  3  4 


TABLE  6-4  ASSSS  WORKING  STORE  MAP 


TABLE  6-5  ASSSS  PROGRAM  COOING  CHART 


TABLE  6-5  (CONTINUED) 


TABLE  6-5  (CONTINUED) 


TABLE  6-5  (CONTINUED) 


MLR  MRR  MULTIPLY  ALIR  ACIR  ARIR  ADD  ALOR  ACOR  AROR  AEOR 


L-A.?.Ei.  He  uftrr  HLR  MRR  MULTIPLY  TaLI R  IaCIR  IaRIR  I  ADD  ALOR  IaCOR  IaROR  j  AEOR 
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CHAPTER  7 

SCALE  AND  ADD  THREE  ARRAYS  (SCLA3)  MACRO 

FUNCTIONAL  DESCRIPTION 

SCLA3  MACRO  scales  an  array  X  by  a  constant  A,  scales  an  array  Y  by  a  constant 
B,  scales  an  array  Z  by  a  constant  C,  and  combines  the  three  scaled  arrays  on  an 
element-by-element  basis  to  form  an  output  array  W.  X,  Y,  Z,  and  W  must  be  of  tr.e 
same  number  of  elements  N. 

SCLA3  is  intended  to  provide  arbitrary  linear  mixing  of  three  inputs,  and  for 
convenience,  the  coefficients  A,  B,  and  C  are  accessed  by  the  AE  from  Working  Store. 
Separate  address  registers  are  used  for  the  four  arrays,  and  the  output  may  over¬ 
write  any  of  the  three  inputs  or  the  scaling  parameters  A,  B,  and  C. 

MATHEMATICAL  DESCRIPTION 

WisA*X.+B*Yi+C*Zi;  i-1.2,  ....N 
where  A,  B,  and  C  are  scalers,  and  X^ ,  Y^ ,  Z^  and  W.  are  the  ith  elements  of  arrays 
X,  Y,  Z  and  W,  respectively. 

IMPLEMENTATION 

SCLA3  is  implemented  with  an  eight-instruction  main  loop  which  processes  one 
pair  of  X's,  one  pair  of  Y's  and  one  pair  of  Z's  to  form  one  pair  of  W's.  A 
twelve-instruction  preamble  loads  the  three  coefficients  from  Working  Store  using 
Address  Register  ARO,  and  loads  the  pipe  with  the  first  pairs  of  inputs  from  the 
three  input  arrays.  In  addition,  the  preamble  decrements  the  output  Address 
Register  AR4  by  one  to  allow  the  first  output  to  be  written  into  the  initial  value 
of  AR4.  Each  pass  through  the  main  loop  then  stores  one  pair  of  W's  in  the  output 
array.  The  main  loop  is  controlled  by  BNZ9. 

SCALING 

See  Table  7-1. 
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ARITHMETIC  ELEMENT  CONTROLLER  (AEC)  REGISTERS 
See  Table  7-2. 

ADDRESS  REGISTERS 


ARO  -  The  increment  code  is  set  to  1.  The 
to  the  location  of  the  parameters  A, 
AR1  -  The  increment  code  is  set  to  1.  The 
to  the  first  address  of  the  X  array. 
AR2  -  The  increment  code  is  set  to  1.  The 
to  the  first  address  of  the  Y  array. 
AR3  -  The  increment  code  is  set  to  1.  The 
to  the  first  address  of  the  Z  array. 
AR4  -  The  increment  code  is  set  to  1.  The 
to  the  first  address  of  the  W  array. 


address  field  is  initialized 
B,  and  C. 

address  field  is  initialized 
address  field  is  initialized 
address  field  is  initialized 
address  field  is  initialized 


INCREMENT  REGISTERS 

Not  applicable.  Wrap  codes  must  be  1024  for  INC  0  to  4. 


BRANCH  REGISTERS 

BNZ9  -  Branch  and  Count  Register  controlling  the  main  loop.  Reset  and 
Count  fields  are  initialized  to  N/2-1 ,  where  N  is  the  number  of 
elements  in  each  of  the  arrays.  Branch  Address  field  is  set  to 
SC3LP. 

SINE/COSINE  DESTINATION  REGISTERS 
Not  applicable. 

DECIMATE  REGISTER 
Not  applicable. 


SCALE  FACTOR  REGISTERS 

SFOL  -  Left,  center  and  right  prescalers  -  True/Direct.  Postscaler  - 
True/Left  4. 


ARITHMETIC  ELEMENT  REGISTER  MAP 
See  Table  7-3. 
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WORKING  STORE  MAP 
See  Table  7-4. 
PROGRAM  CODING  CHART 
See  Table  7-5. 
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SCLA3  WORKING  STORE  MAP 
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SCLA3  PROGRAM  CODING  CHART 
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CHAPTER  8 

FILTER  TWO  AND  SUM  THREE  ARRAYS  (F22S3)  MACRO 

FUNCTIONAL  DESCRIPTION 

F22S3  MACRO  provides  filtering  of  two  input  arrays  X  and  Y  to  produce  filter 
outputs  F  and  G,  respectively,  and  then  combines  F  and  G  with  a  third  array  S  to 
form  a  new  output  array  S.  The  two  filters  for  the  X  and  Y  inputs  are  two-pole 
recursive  with  separately  specifiable  gain  and  pole  position  parameters.  The 
parameters  and  the  initial  states  of  the  two  filters  are  obtained  from  Working 
Store.  The  final  states  of  the  filters  are  returned  to  the  same  locations  in 
Working  Store.  Separate  Address  Registers  are  provided  for  each  of  the  input  \ 

arrays  and  the  output  array.  The  output  array  may  overwrite  any  of  the  input 
arrays  or  the  parameters,  but  not  the  filter  states.  All  arrays  must  be  of  the 
same  number  of  elements  N.  The  output  S  may  be  passed  through  a  final  filter 
consisting  of  a  single  zero  at  Z  =  +1  (i.  e.,  a  zero  at  zero  frequency).  This 
option  is  controlled  by  a  Scale  Factor  Registers  1H  and  1L,  with  a  "0101,  0010" 
giving  no  zero  and  a  "0111,  0111"  enabling  the  filter. 

MATHEMATICAL  DESCRIPTION 

The  two  filtering  operations  are: 

Fi  =  U  *Xi  +  A*  Fi-1  +  B  *  Fi-2’  1  =  ■>’  N 

and 

G1  *  V  *  Y.  *  C  *  G1_1  +  D  *  Gi-2;  i  •  1 . N 

Fq,  Gq,  F_j,  and  G_^  are  initialized  by  the  preamble. 

The  summing  operation  is: 

Si  =  Si  +  Fi  +  Gi 1  i  =  1  >  •  •  • »  N- 

X.,  Y. ,  and  S^  are  the  ith  elements  of  the  X,  Y,  and  S  arrays  respecti vely,  and  F 
and  G^  are  partial  results  resident  only  in  the  AE  registers.  F^_^ ,  GN_^ ,  F^  and 
G^  are  returned  to  Working  Store,  overwriting  the  initial  states. 
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IMPLEMENTATION 

F22S3  is  implemented  with  a  twelve-instruction  main  loop,  a  twenty-two 
instruction  preamble  and  a  six-instruction  clean-up  sequence.  The  preamble  loads 
the  parameters  from  Working  Store  via  Address  Register  ART ,  loads  the  initial 
filter  states  via  Address  Register  AR5,  and  loads  the  pipe  with  the  first  pairs  of 
X's,  Y's  and  S's.  Each  pass  through  the  main  loop  then  produces  a  pair  of  filtered 
F's,  a  pair  of  filtered  G's,  and  combines  these  with  a  pair  of  input  S's  to  produce 
a  pair  of  output  S's.  The  clean-up  sequence  returns  the  final  filter  states  to 
Working  Store  and  stops.  The  main  loop  is  controlled  by  BNZ8. 

SCALING 

See  Table  8-1 . 

The  filter  coefficients  A,  B,  C  and  D  are  treated  as  HI 4  numbers.  The  gain 
coefficients  U  and  V  are  treated  as  HI 2  numbers.  Assuming  the  inputs  X  and  Y  and 
the  filter  states  F  and  G  to  be  HO,  the  partial  results  AF,  BF,  CG  and  DG  then  are 
scaled  FI 4,  and  the  modified  inputs  UX  and  V Y  are  scaled  FI 2.  SFOL  is  used  in 
forming  new  filter  states  and  shifts  the  AF,  BF,  CG  and  DG  terms  right  by  two  on 
input.  Outputs  saved  in  ACOR  and  AROR  do  not  pass  through  the  postscaler,  and, 
therefore,  are  F12.  Outputs  passed  back  to  the  multiplier  are  shifted  left  by  four 
to  regain  a  scaling  of  HO. 

The  summing  operation  uses  SFOH  which  right  shifts  the  S  input  four  to  agree 
with  the  F  and  G  inputs,  and  shifts  the  output  left  four  to  achieve  a  final  output 
scaling  of  HO.  The  filter  states  passed  to  the  output  to  be  saved  for  the  next 
execution  of  F22S3  are  also  shifted  left  four  by  SFOH. 

ARITHMETIC  ELEMENT  CONTROLLER  (AEC)  REGISTERS 

See  Table  8-2. 

ADDRESS  REGISTERS 

ARO  -  The  increment  code  is  set  to  1.  The  address  field  is  initialized 
to  the  beginning  of  the  input  S  array. 

AR1  -  The  increment  code  is  set  to  1.  The  address  field  is  initialized 
to  the  beginning  of  the  parameters  block. 

AR2  -  The  increment  code  is  set  to  1.  The  address  field  is  initialized 
to  the  beginning  of  the  input  X  array. 

AR3  -  The  increment  code  is  set  to  1.  The  address  field  is  initialized 
to  the  beginning  of  the  input  Y  array. 
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AR4  -  The  increment  code  is  set  to  1.  The  address  field  is  initialized 
to  the  beginning  of  the  output  S  array. 

AR5  -  The  increment  code  is  set  to  1.  The  address  field  is  initialized 
to  the  beginning  of  the  filter  states  block. 

INCREMENT  REGISTERS 

INCO  through  5  -  The  wrap  field  should  be  set  to  1024. 

BRANCH  REGISTERS 

BNZ8  -  Branch  and  Count  Register  controlling  the  main  loop.  Count  and 
Reset  fields  are  initialized  to  N/2-1 ,  where  N  is  the  number  of 
elements  in  each  of  the  arrays.  The  Branch  Address  field  is  set 
to  F22LP. 

SINE/COSINE  DESTINATION  REGISTERS 

Not  applicable. 

DECIMATE  REGISTER 

Not  applicable. 

SCALE  FACTOR  REGISTERS 

SFOH  -  Left  prescaler  -  True/Right  4,  Center  and  Right  prescalers  - 
True/Direct,  Postscaler  -  True/Left  4. 

SFOL  -  Left  prescaler  -  True/Direct,  Center  and  Right  prescalers  - 
True/Right  2,  Postscaler  -  True/Left  4. 

SF1H  -  Left  and  center  prescalers  -  True/Direct,  Right  prescaler  - 
True/Right  1  or  Inhibit,  Postscaler  -  True/Left  4. 

SF1L  -  Left  prescaler  -  True/Direct,  Center  prescaler  -  True/Right  1 

or  Inhibit,  Right  prescaler  -  True/Direct,  Postscaler  -  True/Left  4. 

ARITHMETIC  ELEMENT  REGISTER  MAP 
See  Table  8-3. 

WORKING  STORE  MAP 
See  Table  8-4. 

PROGRAM  CODING  CHART 
See  Table  8-5. 


8-3 


GHO,  GLO 


SHO,  SLO 


NSWC  TR  81-313 


TABLE  8-5  F22S3  PROGRAM  COOING  CHART 


F/G  G/t 


AD-A119  Ml 

UNCLASSIFIED 


naval  surface  weapon*  center  silver  spring  md 
ARITHMETIC  PROCESSOR  (AP)  microprograms  FOR  OASS.(U) 
AUG  *1  R  H  DAVIS*  D  M  GREATHOUSE 
NSWC/TR-81-313  SGI-AO-FSOO  008 


TABLE  8-5  (CONTINUED) 
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CHAPTER  9 

DEMODULATED  NOISE  (DEMON)  MACRO 


FUNCTIONAL  DESCRIPTION 

The  DEMON  MACRO  generates  two  channels  of  modulated  broadband  noise  and  adds 
the  two  on  a  sample-by-sample  basis  to  a  composite  channel.  Two  channels  of  white 
noise  are  input  to  the  MACRO  which  filters  each  with  a  two-pole  recursive  filter 
to  provide  some  spectral  shaping  of  the  broadband  noise.  The  modulating  function 
for  each  channel  is  derived  by  truncating  a  sinusoid  with  arbitrary  amplitude,  mean 
and  frequency. 


MATHEMATICAL  DESCRIPTION 


DEMON  computations  are  conveniently  grouped  into  three  phases: 

a)  Generate  Modulating  Functions 

b)  Generate  Filtered  Noises 

c)  Form  Composite  Output  Signal 
Phase  1 

For  the  first  DEMON  channel  (AEO): 

S^.  =  COS(ARO)  ARO  =  ARO  +  INCO 

For  the  first  channel  (AE1): 

S^.  =  C0S(AR2)  AR2  =  AR2  +  INC2 

where  i  is  the  time  index. 

The  sinusoid  is  scaled  by  K 


and  then  offset  by  L,  -  T,  (Level  -  Threshold)  to  form  a  test  parameter  P 

f  h 

The  modulating  function  value  for  the  itn  sample  is  then  given  by 


li‘ 


V’ 

KSli  +  Li 


li 


<  0 


li 


>  0 
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Similarly,  for  the  second  DEMON  channel  (AEO)  and  the  second  channel  (AE1) 


S2i  =  COS(ARl) 
S2j  =  C0S(AR3) 


KS 


2i 


K  *  S 
k2  ^2i 


DM2i  - 


T2; 

KS 


+  I  * 
2i  u2’ 


AR1  =  AR1  +  INC1 
AR3  =  AR3  +  INC3 


p2i  i  0 


P2i  >0 


Phase  2 

For  the  first  channel,  the  input  white  noise  samples  are  passed  through  a 
two-pole  recursive  filter 

+  C  *  Y, 


Yli  =Ai  *Xii  +B*Yi(1_i) 


'1(1-2) 


where  A1  is  a  gain  coefficient,  B  is  equal  to  the  TWOR  filter  coefficient  and  C 
is  equal  to  the  MRSQ  coefficient.  Similarly,  for  the  second  channel 

Y2i  =  a2  *  x2i  +  B  *  Y2(i-1 )  +  C  *  Y2(i-2) ’ 

Phase  3 

The  filtered  noise  channels  are  multiplied  by  the  modulating  functions 
Dli  =DMli  *Yli 
°2i  =  DM2i  *  Y2i 

and  combined  with  the  composite  channel  to  form  the  final  output 

TGi  =  TGi  +  Dli  +  °2i  ’ 

IMPLEMENTATION 


DEMON  is  implemented  with  a  20-instruction  main  •'m»  a  30-instruction  preamble 
and  a  6-instruction  clean-up.  The  preamble  loads  the  ,  liters  and  filter  states 
from  Working  Store  (WS)  via  Address  Registers  AR6  and  AR7,  generates  the  first  two 
points  of  each  periodic  modulating  signal  via  Address  Registers  ARO  through  AR3, 
loads  the  first  pair  of  channel  inputs  via  Address  Register  AR4,  and  loads  the 
first  pair  of  target  data  via  Address  Register  AR5.  Since  the  target  output 
overwrites  the  target  input,  a  dummy  read  of  -1  on  AR5  is  executed.  Each  pass 
through  the  main  loop  produces  two  pairs  of  filtered  Y's  which  are  multiplied  by 
the  corresponding  modified  modulating  signals.  The  products  are  then  added  to  the 
target  signal  to  produce  a  composite  output  target  signal.  The  cleanup  sequence 
returns  the  final  filter  states  to  Working  Store  and  stops.  The  main  loop  is 
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controlled  by  BNZF. 

SCALING 

See  Table  9-1 . 

ARITHMETIC  ELEMENT  CONTROLLER  (AEC)  REGISTERS 
See  Table  9-2. 

ADDRESS  REGISTERS 

ARO  -  Phase  address  for  first  DEMON  channel  (AEO). 

AR1  -  Phase  address  for  second  DEMON  channel  (AEO). 

AR2  -  Phase  address  for  first  DEMON  channel  (AE1). 

AR3  -  Phase  address  for  second  DEMON  channel  (AE1 ) . 

AR4  -  The  increment  code  is  set  to  1.  The  address  field  is  initialized 
to  the  beginning  of  the  input  X  array. 

AR5  -  The  increment  code  is  set  to  1.  The  address  field  is  initialized 
to  the  beginning  of  the  input  TG  array. 

AR6  -  The  increment  code  is  set  to  1.  The  address  field  is  initialized 
to  the  beginning  of  the  parameters  array. 

AR7  -  The  increment  code  is  set  to  1.  The  address  field  is  initialized 
to  the  beginning  of  the  filter  states  array. 

INCO  -  Frequency  increment  for  the  first  DEMON  channel  (AEO). 

INC1  -  Frequency  increment  for  the  second  DEMON  channel  (AEO). 

INC2  -  Frequency  increment  for  the  first  DEMON  channel  (AE1). 

INC3  -  Frequency  increment  for  the  second  DEMON  channel  (AE1). 

INC4,  6,  7  -  Wrap  code  must  be  set  to  1024. 

INC5  -  Wrap  code  -  1024,  increment  field  -  +2. 

BRANCH  REGISTERS 

BNZF  -  Branch  and  Count  register  controlling  the  number  of  executions 
of  the  main  loop.  Count  and  Reset  fields  are  set  to  N/2-1 , 
where  N  is  the  number  of  elements  in  the  target  array.  The  Branch 
Address  is  DMNLP . 

SINE/COSINE  DESTINATION  REGISTERS 

SCD01  -  Destine  Cosine  only  to  AEO. 

SCD23  -  Destine  Cosine  only  to  AE1 . 
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DECIMATE  REGISTER 
Not  applicable. 

SCALE  FACTOR  REGISTERS 

SFOH  -  Left  prescaler  -  True/Right  4,  center  and 
True/Direct,  postscaler  -  True/Left  4. 
SFOL  -  Left  and  center  prescalers  -  True/Right  2 
True/Direct,  postscaler  -  True/Left  4. 
SF1H  -  Left  and  center  prescalers  -  True/Right  2 
True/Direct,  postscaler  -  True/Left  2. 
SF1L  -  All  -  True/Direct. 

ARITHMETIC  ELEMENT  REGISTER  MAP 

See  Table  9-3. 

WORKING  STORE  MAP 

See  Table  9-4. 


right  prescalers 
right  prescaler 
right  prescaler 
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TABLE  9-1  DEMON  AP  MICROPROGRAM  SCALING 
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TABLE  9-2  DEMON  AEC  REGISTER  MAP 


Incr  Code  Address  Register  MS  Wrap  Increment  Register  Memory  Address  Register 


0  3  4 


15  0  3  4 


ANGLE  01 

ANGLE  02 

FREQ  01 

FREQ  02 

Sin/Cos  Destination  Register 
AEO  AE1  AE2  AE3 


ANGLE  11 


ANGLE  12 


1  XI,  X2 


1 


FREQ  11 


FREQ  12 


1 

Parameters 

1 

States 

Decimate  Register 
Count  Reset 


Comments 


BR/BCR  Unconditional/ 

Conditional  Branch  BNZ  Branch  and  Count  Registers 

Branch  Address  Count  Reset  Spare  Branch  Address 

0 _  15  0  78  1516  21  22  31 
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TABLE  9-3  DEMON  AE  REGISTER  MAP 


M 

LR 

MRR 

UH 

K1 

A1 

01  0H 

s 

DM 

1H 

K2 

A2 

1 L  1H 

B 

c 

2H 

Yll 

Y1 2 

2L  2H 

XII 

X 1 2 

3H 

Y21 

Y22 

3L  3H 

X21 

X22 

ALIR 

ACIR 

(JH 

B 

Y1 

0L  0H 

c  1 

Yl" 

1H 

TGI  | 

TG2 

1L  1H 

c  i 

Y2" 

2H 

L1  1 

T1 

2L  2H 

0  1 

3H 

L2  1 

T2 

3L  3H 

DU^MY 

J9L 

1L 

2L 

3L 


91 
1 L 
2L 
3L 


® _ .  _  AEOR 


AX1 

0H 

Yl  1 

I 

Yl  2 

AX2 

1H 

Y21 

1 

Y22 

KS 

2H 

TGI 

[ 

TG2 

D2 

3H 

i 

0L 
1 L 
21 
3L 


0L 

11 

2L 

3L 


ALOE - ,  ,  ACOR  MOP 


Yll 

9 

T1 

9 

niAUP 

p 

Y21 

T2 

1 

9 

Y12 

2 

2 

Y22 

3 

3 

| 
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DEMON  WORKING  STORE  MAP 


DEMON  PROGRAM  CODING  CHART 


TABLE  9-5  (CONTINUED) 
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BBBBBbBHBB 

RB 

BB 

BB 

BP 

88 

BB 

BB 

BB 

BB 

BB 

BB 

BB 

PB88BB88B8 

H88B88H888 

88 

88 

88 

BB 

BB 

BRBBbBBBB 

BBBB888BH 

BB 

BB 

BB 

MB 

BB 

•COMOECK  BBAMP 

• 

• 

• 

eB  6ENEPATI0N  12/05/75  WHO 

18  STEP  PREAMBLE 

BBAMP 

AECW 

LOGxO 

Ml  =  NOT  0 

A 

AECW2 

RDWA*MLR.2.AL1R.?.AR.1 , 
tomlr.oh 

READ  AVI .  COUNT  1 

Ml 

*B 

AECW 

LOGxO 

Ml  =  NOT  0 

0 

AECW2 

RDWRrAL I R . 3 • AW . 1 ,AAR, 
tomrr.oh 

READ  0 >/ 1 

Ml 

*1 

AECW 

SE.OH 

NOP 

2 

AECW2 

RDWBxMWR.O.AR.O, 

MLR.OH.MRR.OH 

READ  A  1 , A? 

PI  x  M 1 *M 1 

•3 

AECW 

ALIR.3.SE.0M 

OAV  *  D V 1 

A 

AECW2 

TOACIR.O.TOAROR.O 

PI.  OAV 

S 

AECW 

MLR.?H,MRR.OH, 

AL/R.2M.AROP.O,SE.OH 

VA 1  x  A V 1  * A  1 

CAV?  *  AVI  .  DAV 

•6 

AECW2 

R0WR=MLR.3,AR.3,SAH, 

AClR.OL.SF.OH, 

toacor.o, tomlr. ol 

BACK  UP  AP3 

PI  *  LOW  PART 

CAV2 

•  7 

• 

AECW 

ALIR.2LC, SE.OH, 

TOAPIR.O, 

TOAROR.l .POUND 

CT  =  -COUNT 

VA  1 

PI,  ROUND  CAV2 

*8 

• 

AECW2 

ROWBsHRR. 1 ,AR.2« 
MLR.OL.MRR.OL, 

TOACOR.l, 

ARIR.O.SF.OL 

HEAD  R 1 

VA?=CAV2.A2 

CT 

PASS  V A 1 

•  9 

• 

• 

AECW 

T0MLP.1H 

VA  1 

10 

AECW2 

P0WB=MWR.2.AR.?,AAR, 
TOAPIR.O. ROUND 

READ  R2 

VA2, HOUND  VA 1 

*11 

AECW 

MLR.1H.MRR.1H, 

ARlP.O.SF.OL 

BHRlxVAl.RRl 

PASS  VA2 

*12 

AECW2 

ROWBxMRR.O.AR.O.AAR, 

mlp.ih.mrr.il. 

ACOP.O.AROR.O.SF.OH, 

TOMLR. il 

RE AO  A3.AA 

BHIlxVAl.PIl 

CAV3=CAV2»DAV 

VA2 

*13 

• 

• 

AECW 

TOALIR.O. 

TOACOR . 0 , TOMLW . OH , ROUND 

6HRI 

CAV3. ROUND  VA? 

*  1 A 

AECW2 

MLR.lL.MRR.2H, 

TOAlIW.1 .BOUND, 
ACOP.O.AROR.O.SF.OH 

BHR2xVA2.RR? 
eHIl, ROUND  CAV3 
CAV4xC4V3*0AV 

*15 

• 

AECW 

MLR.OH.MRR.OH, 

TOACOR. 0. TOMLR. OL 

VAjxCAV3.A3 

CAVA 

•18 

AECW2 

MLR.1L.MRP.2L, 
TOACIR.O.TOAEOR.l. ROUND 

8HI ?x VA2.P I 2 

BHR?,8l «0 .ROUNO  CAV4 

*17 

A-2 
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10  INSTRUCTION  MAIN  LOOR  WITH  FMHEDOtD 
S  INSTRUCTION' VARIATOR  FUNCTION  CHANGE 


AECW 

mL».oL,mpp.0L, 

VAa  =  CAV<*.  Aa 

•18 

TOArip.o. 

VA  3 

• 

ACOR.l.AROR.l, SF.OH, 

CT=CT  *R 1 

• 

TOAEOP. 1L, POUND 

Mil,  ROUND  BP1 

AECW2 

PDWS*MRH.1,A«.2.AAP, 

READ  P3 

*1* 

TOACIR.l, ROUND, 

BHI2, POUND  HU 

• 

ARIP.O.SE.OL, 

PASS  VA3 

• 

TOACOR. 1 ,SETCS 

CT , SET  STATUSft 

AECW 

TOARIR.O, 

VAa 

*20 

AClR.0,ALIR.0t,SE.0L, 

BR?=hHR?-BHR1 

• 

tomlr. ih 

VA  i 

AECW2 

RDWRsMRR.2, AR.2.AAR, 

READ  Ra 

*21 

ARIP.O.SF.OL, 

PASS  VAA 

• 

TO AEOR.2H, ROUND, 

MR?, ROUND  VA3 

• 

BCR.fl 

BRANCH  ON  ST6 

AECW 

MLR.1H,MR«.1H, 

BHR3=VA3.PR3 

*22 

tomlr. 1L. ROUND, 

VAa. ROUND  BR2 

• 

ALIR.I.ACIR.IC.SE.OL 

B I 2aBHI 1 -BHI2 

AECW2 

RDWRaMRR.n, AR.O.AAR, 

READ  AS.A6 

*23 

MLR.lH.MftR.lL. 

BHI3*VA3.»13 

• 

TOAEOR.2L. ROUND, 

MI 2, ROUND  VAA 

• 

ACOP.O.AROR.O. SF.OH 

CAVS*CAVA»DAV 

VARIATOR  FUNCTION  CHANGE 

EXECUTES  ONLY  WHEN  CT  s  0 

AECW 

SF. oh, ROUND 

NOR, ROUND  MI 2 

24 

AECW2 

RDWMsAL Ift.?,AR,l .AAR 

READ  AV, COUNT 

25 

AECW 

sf.oh 

NOP 

26 

AECW2 

ftOWH=ALlH.3,AR.l,AAft 

READ  DV 

27 

AECW 

SF.OH 

NOP 

28 

AECW2 

ALIP.2LC, SF.OH 

CT  =  -COUNT 

29 

AECW 

MLR . 1 H ,MWR . 1 H , 

BHR3=VA3*PP3 

*30 

ALIP.3.SF.0H. 

DAV  =  DV 

• 

toacow.i 

CT 

AECW2 

mlp.ih.mrr.il. 

BHI 3=V A3*R I 3 

*31 

ALIR.2H, SF.OH, 

CAV  *  AV 

• 

toaror.o 

DV 

CONTINUATION  oe  MAIN  LOOP 

AECW 

TOALlR.O, ROUND, 

HHR3.  ROUND  Bl2 

*32 

TOACOR.O, TOMLR. OH 

CAV5 

AECW2 

MLR.iL.MRR.2H, 

bhra=v aa.rra 

*33 

toalir.i, round. 

BHI3, POUND  CAVS 

• 

ACOO.O.AROR.O.SF.OH, 

CAV6=CAVS»0AV 

• 

WRMaAEOO, 1 .AR.3.AAR, 

WRITE  dl 

• 

HNZ.R 

6R  TO  MB1 

AECW 

MLR.0M.MRR.0H, 

VA5=CAV5.A5 

*34 

TOACOR.O, TOMLR. OL, 

CAV6 

• 

ALIR.O.ACIR.OC.SF.OL 

BP3=BHR3-BHR2 

AECW2 

MLR.1L.MRR.2L. 

BHlAaVAA.RlA 

*35 

TOAFOO.IH, ROUND, 

6R3. ROUND  CAVS 

• 

ALIO.lC.ACIR.l.SF.OL, 

HI3=BHI?-BHI3 

• 

TOACIR.O, 

HHR4 

• 

WRMaAEOR,?  »  AR, 3.  A AR 

WRITE  d? 

OUTSIDE 

LOOP  FOR  SIZE  GREATER  THAN  51? 

AECW 

TOAFOR.1L. ROUND 

BI3. ROUND  8R3 

36 

AECW2 

BNZ.9. ROUND 

BP  TO  MB  1, ROUND  BI3 

37 

AECw 

MLR.0H.Mftp.0H 

VAS=CA VS. AS 

38 

AECW2 

MLP, 1L.MRR.2L, 

WHIa=VAA,RI4 

•39 

ALIP.1C.ACIP.1 .SF.OL 

H I  3=eH 1 2-MHl 3 

AECW 

SF.OH 

NUP 

AECW2 

BO. A 

BRANCH  TO  SMU 

AECW 

SF.OH 

NOR 

AECW2 

sf.oh 

NOP  CONTINUE 

A-3 


NSWC  TR  81-313 


000000000 

00000000 

DDODODOD 

00 

00 

DO 

OODOODOODO 

OOOOODOOOO 

OODOODOODO 

00 

000 

DO 

r>n 

00 

00  D 

DO  D 

DO 

00 

0  DO 

DO 

00 

00 

DO 

00 

00 

00  00 

00 

00 

uoooooooo 

DP 

00 

00 

00  OD 

on 

00 

00000D00D 

DO 

DO 

00 

00  00 

00 

00 

DO 

DO 

00 

DO 

00  DO 

00 

00 

0  00 

00  0 

DO 

00 

0  00 

nooouooooo 

OOOOOODOOD 

DODODDODDO 

OODOODOODO 

DO 

000 

OOOOUOOOO 

OUODOOOD 

00000000 

0000000000 

on 

DO 

•coMoecK 

■  OSCLN 

OSCLN 

aecu 

SE.OH 

NOP 

-? 

AECw? 

L0G=0* 

CREATE  MINUS  1 

*-l 

ROWflsAW.?,AR.S, 

read  eama 

• 

SCMR.?HsAP.l 

GET  [51  PARAMETER 

AECW 

L0G=0 , TOMRR .Oh, TOAROR . 0 

OES  MJ 

0 

AECw? 

RDWfisMLP.O.AR.S.AAR, 

READ  es.a 

*1 

TOMER. 1H 

OES  Ml 

AECW 

SE.OH 

NOP 

? 

AECW2 

ROWRsACIR. 1 , AR.fr, 

READ  COOP 

*3 

SCMR.2H3AR . 0 • 

GET  I 7 f  PARAMETER 

• 

«LR. lH.MRR.OH 

P  1  sM  [  *M [ 

AECW 

SE.OH 

NOP 

4 

AECW? 

PDWf?=MLR.3.AR.4,TOALlR.O 

RE AO  SL.P.OES  PI 

s 

AECW 

ACIR.1 .SE.OH 

PASS  CoDP 

6 

AECW? 

ALlR.OL.SE.OH.TOALOR.O, 

P1=L0(P1) ,  DES  COa 

*7 

ROWPsMHW.O.AR.? 

READ  EM 

AECW 

ALIP.OL, SE.OH, TOMLR . 1H. TOACOR , 3 

PlsLO(Pl) ,OES  PI 

AECW? 

T0AL0R.1,T0MRR.3L, 

OES  PI,  PI 

*9 

ACOP . 3 , AROR . 0 , SE . 3L  « 

LMsMl.Pl (R?) 

• 

ROWP=ACIR.O.AR.4,AAR 

READ  COO 

AECw 

mlr.ih,mrr.?h. 

7*7*P1 

•  A 

TOACOR. 3 

DtS  LM 

AECW? 

MLR.iH,MRR.?L, 

5=S*P  I 

ALOR.O, ACOR.3»L0GsN 

OasCOOP  AN0  LM 

AECW 

toacir.2»toaror.o. 

DES  7. OES  D4 

*C 

mlr.om.mrr.oh. 

ESM=ES*EM 

• 

ACIR. 0, SE.OH 

PASS  CDO 

AECW? 

MLR.-JH,MRR.3L,T0ACIR.3, 

SL sPl *SL ,  DCS  5 

*0 

ACIP.?L,SE.?H,T0AL0R.3, 

7»L0(7J*?.0ES  C DO 

• 

ROWR=MRR.l  ,AR.3 

READ  am 

AECW 

TOALJR.?, 

DES  ESM 

•E 

AClR.3L.SE.OH,TOALOR.2 

S=L0<51 .  OES  7 

AECW? 

TOACJR. 2, TOACOR.?, 

DES  SL .  S 

*E 

alip.?,se.ih 

ESMasESM (P4) 

AECW 

TOAPOR.l 

DES  E5M4 

G 

AECW? 

MLP.1L.MRR.3L, 

PPsPvRl 

*H 

ACIR. 0. APOD. l.SE.lh 

NCD4S(COO»ESMA> |R?) 

OSCLP 

AECW 

TOAPOR.?,TOMLR.1L, 

OES  NCD4.NC0 

*10 

ALOR.3,ACIR.2.SE.0h 

CD0=Cl>O*SL 

AECW? 

toaeor.i. 

OES  COO 

*11 

ACIR. JH,AROP.O,SE.3L, 

cnPAsH(cnoP) »da 

• 

ROWR*AR.O.AR.6.TOALIR.3 

REaO  CD.  DES  pp 

AECW 

ACOP, 3 , AROR, 2, LOG* N, 

OD=CD  AND  LM 

*1? 

TOACOR. 0» 

DtS  C0P4 

• 

A-A 


"L 
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AECw? 

MLH.1L.Mhw.3L 

NC*P)*NCD 

AlIR.3.ACoR.OC,AHOH.OC,SF.OL, 

TOAfOH.?, 

WHMsAEOH. 1 ,AR.4 

CCsPP-CDR4-D4 

DES  00 

WRITE  CDO 

*13 

• 

AECw 

AROp. 1 .ACIW.O.SF .OH, 

toalih.o.toacuh.i 

C=CD0*FSM4 

DES  NC,  CC 

*14 

AECW2 

ALlP.OL.ACOR.l .AROR.2.SF.3H, 
R0wq=AR.2,AR.4,SAR, 

SCMr.oHsAH.1 ,TOAEQR.?h 

R=L(NC)/?*CC»NCD4 
READ  SL,  P 

COS.SlN  (0) ,  DES 

•IS 

• 

C 

AECw 

ACIR.2L.SF .0H.T0AE0R.3 

SL=LO(SL) ,  DES  P 

lfc 

AECW2 

ALOR.2,SF.fiH,TOAEOH.3H, 

wRM*a£0H.?,AR.6 

PASS  7,  DES  SL 
WRITE  CD 

*17 

AECW 

MLR.0C.MRR.1M, 

AC0R.2.SF.0H, T0MLR.3H 

A*A*AM 

PASS  S.  DES  7 

*18 

AECW? 

T0MLR.3L. 

WRMsAE0H.3,AR.4. 

SCMp. im=ar. 3 

DES  5 

WRITE  SL,  P 
COS.SlN  (P) 

*19 

• 

AECW 

MLR.3H.MRR.0H, TOACIR.O 

7*C»  DES  A 

1  A 

AECW? 

MLR. 1H.MRR.0H. 

AClP. 0  «  SF  «1C« 

ROWRsAL IR.O.AR.O.AIRO 

P1*C=C 

PASS  A 

RE  AO  BM2 

*1B 

• 

AECW 

MLR. 1H.MRR.0L, TOACIR.O, 

TOMlb.OL 

S=P1*S,  DES  7C 

DES  A 

*1C 

AECW? 

MLR. 3L ,MWh. OL.TOACIR.1. 
ACIO.OA.SF. OH, WOUND, 

ROWRsAL 1 R, 1 , AR.O , AAR 

5*S.  DES  C 

PASS  7C ,  RND  A 
READ  HM 1 

•ID 

• 

AECW 

MLR. OL.MRR.IH, TOACIR.O, 

AC  Ip . 1L A , SF . OH , TOMLR . 2H 

A*K .  DES  S 

C«LO(C) ,  DES  7C 

*1E 

AECW? 

AC Ip. OL A. SF. OH, ROUND. 

TOMLR. 3H.T0AC1R.1, 

ROWRa AL IR. 2, AR.O, AAR 

S=LO(S) ,  RND  7C 
DES  C.  5S 

READ  BO 

•IF 

• 

AECW 

MLR, oL.MHr.IL, TOACIR.O, 

AC Ip. 1A.SF.0H, TOMLR. 3L 

A»2.  DES  AK 

PASS  SS,  OES  S 

*20 

AECW? 

AC  I p . 0 , SF . 2L , T  OMLR . 2L . 

ROWRs AL IR. 3. AR.O, AAR 

PASS  AK,  DES  5S 
REAU  BP1 

*?\ 

AECW 

MLR. 3H, MRP. OL. TOACIR.O, 

T  OACOR . 0 , TOMRr . 2h , ROUND 

C*S,  DES  A2 

DES  AK.  RND  SS 

*?? 

AECW2 

MLR . ?L , MWH . OH , ROUND , 

ACIp.0,SF.?L, 

RDWRsAR,2,AR.5,AlRt 

5S*C  .  ROUNO  AK 
PASS  A J 

READ  Fma,  AMA 

*?3 

• 

AECW 

MLR. PH.MHR.2H, TOACIR.O, 
TOACOH.1. TOMRR. ?L 

7C*AK,  DES  CS 

DES  AZ 

*?4 

AECW? 

MLR. 3L,MRr. OL.TOACIR.1  , 

AC lR.0A,SF.?L, ROUND, 
R0WR=MLH.0,AR.5,AAW 

S*5,  DES  5SC 

PASS  CS,  RND  A? 
READ  FS»  A 

*25 

• 

AECW 

MLR. ?H,MRR.2L, TOACIR.O, 

AC IR . 1 . SF . ?H , TOMRR . 3H 

7C*AZ.  DES  7CAK 
PASS  SSC,  OES  CS 

*20 

AECW  2 

MLR.3L,MRR.?H,TOACIP.I, 

AC I P . 0 . SF . ?L , ROUNO  « 

TOMlp. 1L. 

RDWRaMHR.O, AR.2 

AK*S,  DCS  SS 

PASS  7CAK ,  RNO  CS 
DES  5SC 

READ  Fm 

♦27 

• 

• 

AECW 

MLR. 3L.MRH.2L, TOACIR.O, 

ACIP. 1  A, SF.2L, TOMLR. ?H, ROUND 

A2*S,  OES  7CA 2 

PASS  SS.DES  7CAK.R 

*28 

SSC 

AECW? 

MLR.3H,MRR.2H,TOARlR.o, 
AClP.O.SF ,?L,T UMRW . 3L . ROUND . 
ROWPsMHH. 1 , AH. 3 

AK*C,  DES  ASK 

PAS  7CAZ.DES  SS,R 
HEAD  AM 

*29 

7CAK* 

AECW 

MLR.?H,mHh.3H,T0ARIR.1 , 

ACOP.OC, AHIR.o.SF ,2H, 

7CAK*CS.  DES  AS 7 

T  1  = ( ASK-AK ) /2 

*2  A 

• 

f 
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Kj 

TOMt_P  ,  2L »  ROUND 

UES  7CAZ,  RNO  SS 

S  AECi»2 

MLR.3H.MRR.2L.T0ARIR.?, 

AZ*C,  OES  AC* 

•?a 

■ 

ACOR. 1C. ARlR. 1 .SF .2H, 

12= 1ASZ-AZ) /<- 

• 

■ 

tqapor.o. round. 

UES  Tl,  RNO  7CAZ 

• 

■  ' 

R0wr=MLR.3,AR.4,AIRE 

Rt  AO  SL,  P 

1  :  AECW 

MLR.2L.MRR.3H.TOACIR.0. 

7CAZ*CS.0ES  7CCSA* 

*2C 

ft  , 

ACOP.0.ARIR.2C.SF ,?h. 

73*  ( AK-ACM  /2 

« 

toaror.i 

DES  T? 

B  AECW2 

Al»  I R.  OH,  Ac  1R.0C.SE.  OH. 

h?R=BM?R-7CCSAK 

*20 

■ 

T  OAR I R . 3 . T  OAROR . ?  . 

DES  ACZ.  T 3 

• 

ft 

ROWR=ACIW.O.AH.4,AAR 

RE  AO  COO 

1  AECW 

MLR.  IU.MRR.2H,  TC  *ClR.  1  , 

5SC*AK ,0E S  7CCSAZ 

*2E 

ft 

ACOR . 1 , AR I R , 3C , SF ,2H, 

T4  = ( AZ-ACZ ) /2 

• 

ft 

toaeow.oh 

OES  M2S 

ft  AECW2 

al  ir.ol.acir.ic.se. OH, 

M2l=BM?l-7CCSAZ 

*?E 

ft 

T OAPOR • 3  *  ROUND , 

DES  74, RND  M2R 

• 

ft 

RDWRaACIR. 1 , AK.6,  AAR 

READ  CO 

I  AECW 

MLR.1L.MRP.2L.T0ACIR.2, 

SSC*AZ,  OES  SSCAK 

*30 

ft 

ACIP.O.SE.OH.TOAEOR.OL 

PASS  COO.  DES  M? 1 

I  AECW2 

ALIR.1H,ACIR.2C,AROR.O, 

M1R=HM1R-5SCAa*T1 

*31 

ft 

SF.?h.TUAL0R.3.R0UND, 

DES  CDO,  RNO  M2I 

• 

ft 

Rewire 1R.0.AR.0.AAR 

READ  BP2 

1 

I  AECW 

ACIR.1.SE.OH,ToACIR.3.TOA£OR.1H 

PASS  CO, DES  SSCAZ.M1R 

32 

I  AECW2 

ALIR.1L,ACIR.3C,AR0P.I.SE.2h, 

Ml  I=BM1 1-5SCAZ.T? 

*33 

1 

I 

TOALOR.O, round. 

DES  CD4.RND  MIR 

» 

1 

ROWB*ALIR.l .AH.©, AAR 

READ  BP3 

| 

1  ,  AECW 

ALIR.2H,ACOR.O,AR1R.?,SE.?H, 

0R=R0R*AK/2*ACK/2 

*34 

j 

1  1 

T0AE0R.1L 

OES  Mil 

1  AECW2 

ALIR.2L.AC0R. 1.APIR. 3, SE.2H, 

0l=8OI»AZ/2*ACZ/2 

*35 

| 

T0AE0R.2H, ROUND, 

DES  OR.  RND  mu 

• 

I 

WRMsaEOR.O.AR.O.AIWE 

WHITE  M2 

1  AECW 

ALlR.3H,ACOR.OC,ARlO.OC.SE.?H, 

PlR=BPlR-AK/2-A5K/2 

*36 

1 

T0AP0R.2L.R0UN0 

OES  01,  RNO  OR 

1  AECW2 

ALlR. 3L» ACOR. 1C.AR1R.1C.SE.2H. 

PI  I  =  BP1 I-AZ/2-ASZ/2 

*37 

T0AE0R.3H, ROUND, 

DES  PlW,  RNO  01 

• 

WRMsAEOR .1 .AR.O.AAR 

WRITE  Ml 

AECW 

MLR.2H.MRR.3U, 

7C AK*5S 

•38 

: ! 

ALlR .OH. ACIR.2 , AROR .? .SE ,?H, 

P2R=BP2R.5SCAK.T3 

• 

TO AE0R.3L, ROUND 

OES  PI  I,  RNO  PIP 

AECW2 

MLR.2L.MRR.3L, 

7CAZ*SS 

*3R 

F 

AL IP. OL , ACIR. 3 t AROR. 3 ♦ SE.2H, 

P2I=PP?I»5SCAZ.T4 

« 

TOAEOR. OH, ROUND, 

OES  P?R «  RNO  P  V I 

• 

WRMsAEOR.?, AR.O.AAR 

WRITE  0 

AECW 

ALOP.l.SF , OH , TOAC IR ,2  * 

PASS  PI,  DES  75SCAK 

*3A 

MLR.OH.MRR.OH, 

ESM=ES*EM 

• 

TOAEOM.OL, ROUND 

DES  P?I.  RNO  P2R 

aecw? 

ALIR.1H.ACIR.2.SE.0H, 

P3R=BP3R*7SSCAK 

•35 

A 

TOAC I R . 3 , TOMRR . 3L , ROUND , 

DES  7SSCAZ.P1 .RND  P2I 

• 

.1 

WRMsaEOR. 3, AR.O.AAR 

WRITE  PI 

AECW 

AHP.1L.ACIR.3.SE.0H, 

P3t=8P3I*7SSCAZ 

*3C 

TOALIR.2.TOAEOR.1H 

DES  ESM.P3R 

j  AECW2 

MLR.3H.mrr.3L , 

SL=P1«SL 

*30 

toaeor.il.rouno. 

DES  P3I,  RNO  P3R 

• 

WRMsAEOR.O. ar.o.aar. 

WRITE  P2 

• 

ALIR.2.SE.1H, 

FSM4*ESM<HM 

• 

bnz.a 

RR  ON  CT  TO  OSCLP 

j  AECW 

ALOR . 0  » ACOR , 3,L00*N, 

04=COOP  AND  LM 

•3E 

toaror.i. ROUND 

DES  ESM4,  RNO  P3I 

A  AECW? 

TOACTR.2.TOARUR.0, 

DES  SL,  04 

*3E 

f 

AClR.O.AROR.l.SF.lH, 

NCD4*CD0*ESM4 

♦ 

* 

MLR.3L.MRR.3L. 

PP*R 1 *P 

• 

WRMsAEOR. 1 .ar.o.aar 

WRITE  P 3 

AECW 

SE.OH 

NOP 

40 

i  AECW2 

i 

STOP 

41 

l 

’ 

i! 

A-6 

i 

^  .  ..  ,t 

■  ‘  'i'  ■ 

w— 
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s  SS 

SS 

SS 

SS 

SS 

SS 

SSSSSSSSSS 

SS 

SS 

SS 

SS 

SS 
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SS 

SS 

SS 

SS 

SS 

srfft 

AECW 

LOG=0 

CREATE  ALL  ONES 

0 

A£CV»? 

TOMLR. 2H. 

MlsMlMJS  1 

*1 

R0WP=AL1R. 1 ,AP.2 

HEAD  GO 

AECW 

L0G=0 

NEED  TWO  MUS 

? 

AECW2 

TOMpp.  1  ri  • 

Ml 

*3 

RDWA*MLH.0.ALIR.2.AR.2. A AW. 

READ  G 1 

• 

SCMp.O=AP.OlNC 

GET  COS, SIN  (INC) 

AECW 

ALlO.lH.SF  .OH 

FRO=GRO 

4 

AECW2 

ALlR. 1H.SF .OH. 

FI0=GW0 

*5 

ROWRrACIR.O.AP.3. 

READ  GN1 

• 

TOAFOR.OH 

FRO 

AECW 

TOAFOP.OL 

FIO 

6 

AECW2 

rowr=mlw.3,ar.s,aaw 

Dummy  READ  FOR  Fn 

7 

AFCW 

ALIp.2L.ACIR.0L,SF.0L 

CIl=(GINl,GIl)/2 

B 

AECW2 

T0MLP.1L, 

CIl 

*9 

ALlP.2H,ACIft.OH,SF.OC. 

AP I  * ( GRN 1 »GR 1 ) /2 

• 

PDwA=ML«.0,AHR.O,AR.2,AAR 

READ  G? 

AECW 

T  OACOR . 0 , ROUND . 

API.  ROUND  CIl 

•  A 

alir.?hc.acir.oh.sf.ol 

.  CR1  =  ( GRN1 -GR 1 )  /  2 

AECW? 

TOMLR. 1H, 

CHI 

*B 

MLR.IL.MPP.OH, 

CCI1=C*CI1 

• 

ALlP.?LC,ACIR.OL,SF.OL, 

AI 1  =  (GIN1-GI 1 ) /? 

• 

POHA=MRH.2.ACIW.O,AP.3.SAR 

REA!)  GN2 

AECW 

TOACOW.l .ROUND, 

All,  ROUND  CR 1 

*C 

MLR.1L.MRR.0L 

SCI1=S*CI1 

AECW? 

TOAPJP. J, 

CCR 1 

*0 

MLR. 1H.MHP.0L, 

SCR  1 =S*CR 1 

• 

WRMsAEOH.O.AR.A. 

wHITt  FO 

• 

SCMP, 0«AP. 0  INC 

COS, SIN  (2*1 NC ) 

AECW 

T0ALIR.3, 

SCI  1 

*E 

MLR. 1H, MRP. OH, 

CCR 1 =C»CR 1 

• 

ALIR.OHC.ACIR.Oh.SF.OL 

CR?= ( GRN2-GR2 ) /? 

AECW? 

TOAl  I  R.2 ,  T  C'MLH.  1H  , 

SCR1 ,CH2 

*F 

ALIp.OL.ACIR.OL.SF.OL, 

CI?=(GIN2*GI2)/’? 

• 

MLR. OH.MhH,  1H 

MGR?=M1*GR? 

AECW 

T0APIR.2, TOMLR. 1L ,ROUNO, 

CCR.CI.ROUNO  CP 

*10 

MLR.2H,MRR.2H, 

MGRN=M1*GRN 

• 

ALIR.?.aRIR.3,ACOR.O,SF.1H 

FR=SCR,CCI*AH 

AECW2 

TOAlIR. 1 , TOAEOR. OH, ROUND. 

MGR,  FR,  POUND  Cl 

*11 

MLR.?H,MRR.2L, 

MG  1 NsM 1 *G I N 

• 

ALIR.3,AR1R.?C.AC0R.1.SF. 1H. 

F  I  =  SC I -CCR* A  1 

• 

RDWA=MLR.0,ALIR.0.AR,2, AAR 

READ  G 

AECW 

TOACIR. 1 .TOAtOR.OL, ROUND, 

MORN ,  FI,  ROUND  FR 

*1? 

MLR.OL.MPP. 1H, 

MGI=M)*SI 

• 

A-7 
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AL1P.2C,AH1R.3L,ACOR.O,SF.1H 

FWN=-SCR-CC1 *AR 

AECW2 

T0ACI«.2. TOAEOW. 1H, POUND. 

MO IN,  FRN,  POUNO  FI 

*13 

MLR. 1L.MRW.0H, 

CCI=C*CI 

• 

ALIR.3.aPIP.?L,ACOR.1C,SF.1h. 

FlN=SCl-CCW-Al 

• 

ROWAxMWR.2.ACIR.O,A9.3.SAW 

READ  ON 

AECw 

TOAL IR. 2. TOAEOW. 1L. ROUND, 

MOI .FIN, POUNO  FRN 

*U 

MLR.1L.MHW.0U. 

SCl=S*CI 

• 

ALIR. ILC.ACIR.1LC.SF.0L 

AH= ( -MGH-MGRN) /? 

AECW? 

TO Ap I R . 3 . TO4C0R . 0 . HOUND . 

CCI , AR, ROUND  F IN 

*15 

MLR.1H.MHW.0L. 

scw=s*cw 

• 

ALI9.2L.ACIR.2LC.SF.0L. 

A I s ( MG  I »MG IN ) /? 

• 

WRM2AEOR.O.AP.4.AAR, 

WRITE  F 

• 

SCMp.0=AW.0lNC,bN2.e 

COS, SIN, LOOP  ON  «N2b 

AECw 

TOA(_IW,3.  TOACOH .  1  . 

SCI .AI 

*1<S 

MLR. IH. MRP. OH, 

CCRsC*CH 

• 

ALlP.OHC.ACIH.OH.SF.OL 

CR= ( GRN-GP ) /Z 

AECw2 

TOALIR.2.TOMLH.1H. 

scw.cw 

•IT 

MLR.0H.MHR.1H, 

MGR=M1 *GR 

• 

ALlP.OL.ACIR.OL.SF.OL, 

CI=(6lN*GI)/2 

• 

WRMxAEOW. 1 .AR.S.SAW 

WRITE  Fn 

AFCW 

SF.OH 

NOP 

18 

AECW? 

BNZ.C 

OUTSIDE  LOOP 

19 

AFCW 

MLR. 1H.MWR.0H 

CCW=C*CR 

1  A 

AECW2 

MLP.nH.MWW.lH, 

MGR=M1*GR 

*18 

ALlR.OL.ACIR.OL.SF.OL 

CI=(GIN.GI 1/2 

AECW 

SF.OH 

NOP 

1C 

AECW2 

STOP 

STOP 

ID 
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•COMDECK  SCLA2 

SCLA2 

AECW 

SF.OH 

NOR 

0 

AECW? 

rowb=mlh.c«*R.2« 

HEAL)  1ST  PAIR  OF 

XIS 

*1 

SCMR.OHsAR.O 

GENERATE  PARAMETER  A 

AECW 

SF.OH 

NOP 

2 

AECW? 

ROWB=MLH.l .AR.3. 

HEAD  1ST  PAIR  OF 

Y1S 

*3 

SCMR.OL=AR.l 

GENERATE  PARAMETER  B 

AECw 

SF.OH 

NOP 

AEC-? 

SF.OH 

NOP 

AECW 

mlr.oh.mwr.oh 

AXH=A*XH 

4 

AECW? 

MLR.9L.MPR.OH, 

AXL=A*XL 

*5 

RDWB=MLR.?.AR.4.SAR 

BACK  DP  7  POINTER 

AECW 

TOALIR.O, 

AXH 

*6 

MLR. 1H.MRH.0L 

HYH=6*YH 

AECW2 

TOALIR.l. 

AXL 

*7 

MLR.1L.MRP.0L, 

BYL=B*YL 

• 

rdwb*mlr . 0 . ar . 2 , A  ar 

READ  NEXT  PAIR  OF 

xts 

AECW 

TOACIR.O 

byh 

B 

AECW? 

RL)WP=MLR.1.AH.3,AAR, 

HEAD  NEXT  PAIR  OF 

Y  [  X 

*9 

TOACIR.l 

BYL 

AECW 

MLR.OH.MHH.0H, 

AXH=A*XH 

*A 

ALIR.O.ACIR.O, SF.OH 

ZHsAXH.BYH 

AECW2 

MLR.OL.MHH.OH, 

AXL=A*XL 

•  B 

ALlR.l.ACIR.l. SF.OH, 

ZL=AXL  *BYL 

• 

TOAEOR.OH 

ZH 

SC2LP 

AECW 

MLR.1H.MRR.0L, 

8YH=B*YH 

*C 

TOAEOR.OL, TOALIR.O 

ZL,  AXH 

AECW? 

ROWfl=MLR,O.AR.2,AAH, 

RE  AO  NEXT  PAIR  OF 

XIS 

*D 

MLR.lL.MRR.OL, 

BYL=Ft*YL 

• 

TOALIR.l 

AXL 

AECW 

TOACIR.O 

BYH 

E 

AECW? 

ROWRsMLR. 1 ,AR.3, AAR, 

REAL)  NEXT  PAIR  OF 

Y(S 

*F 

TOACIR. 1 , BNZ  .8 

HYL ,  BRANCH  TO  SC2LP 

AECW 

MLR.OH.MHR.OH, 

AXH=A*XH 

*10 

ALlR.OH.AClR.OH.SF.OH 

Zh=axh*byh 

AECW2 

WRMsAEOR.O,AR.4,AAH, 

WRITE  A  PAIR  OF  ZtS 

*11 

mlr.ol.mhr.oh. 

AXL=A*XL 

• 

ALlR.l.ACIR.l, SF.OH, 

ZL=AXL*SYL 

• 

TOAEOR.OH 

ZH 

AECW 

SF.OH 

NOP 

1? 

AECW? 

BNZ. 9 

OOTSIOE  LOOP 

13 

AECW 

MLR.0H.MRR.0H 

A  XH= A*XH 

14 

AECW? 

MLR.OL.MRR.OH, 

AXL* A*XL 

•IS 

ALlR.l.ACIR.l, SF.OH 

ZL=AXL*BYL 

AECW 

SF.OH 

NOP 

16 

AECW? 

STOP 

STOP 

17 
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•COMOECK  ASSSS 

*  MLPTH  -  MULTIPATH  RESAMPLE, SCALE  AND  ADD 

• 

*  IN  PHASE  START  SEQUENCE 


WHO  01/26/76 


MLPI 


AECW 

SF.OH 

NOP 

0 

AECW2 

ROWAaMRR. 1 .ACIR.l ,AR.O 

READ  A,  -1 

1 

AECW 

SF.OH 

NOP 

2 

AECW2 

RDWPaALlR.?, AR. 1 

READ  F 

3 

AECW 

TOALOR.O. SF.OH 

0  a-INHIBIT 

4 

AECW? 

RDWAsMLR. 0.ALIR.0.AR.2, 

READ  DHDL  <11 

•S 

TOAROR.O, SF.OH 

0  TO  OC 

AECW 

TOAPOR.3, 

0  TO  1C 

*6 

ALOP.OC, AROR.OC.ACIR. 1LC* SF.OH 

0NE=-1  -EPS 

AECW2 

rowasmlr.o,alir.o.ap.?,aar. 

READ  DHDL  <2> 

*7 

TOAPOR.l » 

ONE  TO  OC 

• 

ALOR.0C, AROR.OC.ACIR. 1LC, SF.OH 

0NE»-1  -EPS 

AECW 

mlr.ol.mrr.ih. 

AL  *  A*OL 

*8 

T0AP0R.2 

ONE  TO  1C 

AECW2 

R0WHaALIH.3,AR.l ,AAR 

READ  0 

9 

AECW 

MLR.0H.MRR.1H, 

AH  *  A*DH 

*10 

T0AP1P.3, 

AL 

• 

ALIR.2C.AR0R.1. SF.OH 

F  3  =  ONE-F 

AECW? 

ARIR.3.SF.0L. 

PASS  AL 

•11 

T0AC0R.2, TOMLR. 1L 

F3 

AECW 

TOARIR.2.TOMRR.OL, 

AH,  AL 

*1? 

ALIR.2.SF.0H 

FI  a  F 

AECW2 

ROWRaALlR. 1 ,AR.3, 

READ  SHSL(l) 

*13 

ARIR.2  »SF . OL .ROUND, 

PASS  AH, ROUND  AL 

• 

TOACOR.O.TOMLR.1H.SETCS 

FI. SET  STATUS  6 

AECW 

MLR.0L.MRR.1H, 

AL=A*DL 

*14 

ALIR.3C,ACOR.OC«AROR.2TS. SF.OH, 

FAa-Fl-0*lC 

• 

tomrr.oh 

AM 

AECW? 

RDWAsMLR .0 . ALIR.0,AR.2»AAR. 

READ  DHOL (3) 

*15 

MLR. 1L.MRR.0L .ROUND, 

3AL*F3*AL, ROUND  AH 

• 

T  OACOP • 3 , TOMLR • 2L 

FA 

AECW 

MLR. 1H, MRP. OH, 

1 AHsF\*AH 

*16 

T0ARIP.3, 

AL 

• 

ALIR.3,ACOR.O,SF.OH 

F2=F i*D 

AECW2 

R0WBaMLR.3,AR,4,SAR, 

DUMMY  READ 

*17 

TOAPIR.O, 

3AL 

• 

ARlR.3*SF.0L»StTCS» 

PASS  AL.SET  status 

• 

T0AC0R.I«T0MLW.2H, 

F2 

• 

bnz.a 

BRANCH  TO  MIPLO 

AECW 

MLR.0H.MRR.1H, 

AHaA«DH 

•18 

TOACIR.O.TOMRR.OL. 

1  AH, AL 

• 
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MIPLP 


ALIB.3C.AC0R. 1C. AROR. 2TS.SE. OH 

E  3=-E2-0* 1 C 

AECW2 

MLR, ?L, MRR. OH, 

4AH=E  4* AH 

*19 

ALlp.lH.ACIW.O . ARIP. 0 • SE . 1 H, 

OHsSH. 1 AH.3AL 

• 

TOACOR. 2. T0MLR.1L, ROUND, 

E 3, ROUND  AL 

• 

BCR.l 

BRANCH  T0I002  ON  E2 

out  phase  start  sequence 

AECW 

SE  .OH 

NOP 

20 

AECW2 

ROW  a*MRR • 1 . AC iH.l.AR.O, 

READ  A,  -1 

*21 

TOALOR.O.SE.OH 

0= INHIBIT 

AECW 

SE.OH.TOAROR.O 

OTO  OC 

22 

AECW2 

ROWAsMLR.O.AL1W.O.AR.2« 

READ  OHOL(l) 

*23 

TOAroR.3 

0  TO  1C 

AECW 

ALOr.OC.APOH. OC.ACIR.ILC.se. oh 

0NE=-1  -EPS 

24 

AECW2 

RDWRsAL IR . 2 , AR . 1 . 

READ  E 

*25 

ALOr.OC.AROR. 0C.ACI9.1LC.se. oh. 

ONE*— 1  -EPS 

• 

TOAPOR.l 

ONE  TO  OC 

AECW 

MLR.oM.MRR. im. 

ah*a*oh 

•26 

TOAboP.2 

ONE  TO  1C 

AECW2 

R0WR=ALIR.3«AR.1.AAR 

READ  0 

27 

AECW 

MLR.oL.MRP.1H, 

AL*A«DL 

•  28 

TOAriR.2. 

AH 

• 

AL1R.2C.AROR.1.SE.OH 

E4=0NE-E 

AECW2 

bdwa=mlr.o.ahp.o,ar.?,aar. 

READ  DHDL (2) 

*29 

ARIR.2.SE.0L, 

PASS  AH 

• 

TOACOR.3.TOHLR.2L 

E4 

AECW 

TOAPI R. 3, TOMRW . OH . 

al.ah 

*30 

ALIO. 2. SE. OH 

F2=E 

AECW2 

RDWfl*ACIR.3.AW.4,SAR, 

DUMMY  read 

*31 

ARIR.3.SE . OL , ROUNO . 

PASS  AL, ROUND  AH 

• 

TOACOR. 1 .T0MLR.2H.SETCS 

E2 

AECW 

MLR.OH.MRR, 1H, 

ah=a*dh 

*3? 

ALIR.3C.ACOR.1C,AROR.2,SF.OH. 

E  3=-E  2-D.ONE 

• 

TOMPR.OL 

AL 

AECW2 

ROWRsAL IR. 1 . AR.3. 

READ  SHSL(l) 

*31 

MLR  .  J>L  .  MRR .  OH  .ROUND  . 

4AHSF4*AH, ROUNO  AL 

• 

TOACOR. 2. T0MLR.1L 

E3 

AECW 

MLR.PH.MBR.OL, 

2AL*E  2*AL 

*34 

TOARIR.2, 

AH 

• 

ALIR.3.ACOR.1.SE.OH 

E 1*E2*0 

AECW2 

ROWAsMLR.O.ALIR.O«AR.2.AAR, 

READ  OHQL ( 3 ) 

*35 

T0APIR.1. 

4AM 

* 

ARIR.2.SE.0L, 

PASS  AH 

• 

TOACOR. 0.T0MLR.1H.SETCS. 

El 

• 

9N2.0 

BRANCH  TO  MOPLQ 

AECW 

MLR.0L.MRR.1H, 

AL*A*DL 

*36 

TOACIR.l. 

2AL 

• 

ALIR.3C.AC0R.0C, AR0R.2TS.SE. OH, 

E4S-0-E 1 ♦ 1 c 

• 

tomrr.oh 

AH 

AECW2 

MLR.1L.MRR.0L, 

3AL*E3*AL 

*37 

ALlR.lH.ACIR.l .ARIR.1 .SE.IH, 

OH*  SH ♦ 2  AL ♦ 4  AH 

• 

T0MLR.2L. TOACOR. 3. ROUNO, 

E4, ROUNO 

• 

BCR.  4 

BRANCH  TO  0101  ON  El 

INPHASE 

PROCESSING  LOOP 

AECW 

MLR, OH. MRR, 1H, 

AH=A*DH 

•  38 

toacir.o. 

1  AH 

• 

A-1 1 
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mIPLO 


MIPL? 


MJPL1 


STOP 

MQPL.P 


MOPLO 


MOPL? 


AECW? 

AEC  W 

AECW? 

AECW 

AECW? 

AECW 

AECW? 

AECW 

AECW? 

AECW 

AECW? 


ALIP.3C.ACOR. 1C,APOR.?TS.SF.OH, 

TOMpp.OL 

MLR. ?L. MUR. OH, 

ALIp.IH.aCIR.O, ARIR.O.SF. 1H. 
T0AC0R.2.T0MLR.  1L, WOUND, 
WRMxAEOR.O, AR.4.AAR, 

SCR.  I 

MLR.?H,MH«.OL, 

T0APIR.2. 

ALlP. 3, ACOR .l.AROR.OTS.SF ,0H« 
TOAfOR.OH 

PDWAxMLR.O.ALIR.O. AP.?, AAP, 
TOAQIR.l, ROUND, 

ARIP.2  »SF .OL » 

toacor.o.tomlr.ih.setcs 

MLR.0L.MRR.1H, 

TOACIR.l, 

ALlP.3C,ACOR.OC,AROR.?TS,SF.OH, 

TOMPR.OH 

POWPaALlP. 1 .AR.3.AAR, 

MLR . J  L . MRR , OL , ROUNO . 

ALIfl.  1L.ACIR.1  .APIR.USF.1H, 
TOACOR .3 . TOMLR ,?L , 

BNZ.8 

MLR. 1M.MRR.0M, 

T0APIR.3, 

ALIP.3,ACOR.O.AROR.OTS.SF.OH, 

TOAFOR.OL 

TOAPIR.O. 

ARIR.3.SF, OL .ROUNO , 

TOACOR. 1 . TOMLR. 2M.SETCS. 

BCR,  0 
SF.nH 

WRMsAEOR.0.AR.4,AAR, 

8NZ.C 
SF  .OH 
STOP 


F3a-0-F2»lC 

AL 

4AHsF  4*AH 
OHaSH,) AH*3AL 
F 3. ROUNO  AL 
WRITE  OUTPUT  OHOL 
BRANCH  ON  F?  TO  1002 
?AL=F?*AL 

AH 

FlaQ.F2*0C 

OH 

READ  OHDL 
4 AH, ROUND  OH 
PASS  AH 
FI 

AL=A*OL 

?AL 

F4=-D-F1 »1C 
AH 

READ  SMSL 

3AL*FJ*AL. ROUND  AH 
0LaSL,?AL*4AH 
F  4 

BRANCH  ON  CT  TO  MIPLP 
1  AHxF 1 *AH 
AL 

F  2eD*F 1 *0C 

OL 

3AL 

PASS  AL, ROUND  OL 
F? 

BRANCH  ON  Fl  TO  1001 
NOP 

WRITE  LAST  OUTPUT 
BRANCH  TO  STOP 
NOP 
STOP 


39 


•40 


41 


4? 


43 


44 


45 


46 

*47 

48 

49 


out  phase  processing  loop 

AECW  MLR.0L.MRP.1H, 

TOACIR.l, 

ALlR,3C,ACOR.OC,AROR.?TS,SF.OH. 

TOMpp.OH 

AECW2  MLR, lL. MRR. OL, ROUNO, 

ALlP.lH,ACIR.l,ARIP.l,SF.lH, 
TOACOR. 3, TOMLR. ?L, 
WRMsAEOR.O,AP.4,AAR, 

BCR. 4 

AECW  MLR. 1H, MRR, OH, 

T0APIR.3. 

AL I R . 3 , ACOR . 0 , AROH .O.S.SF.OH, 
TOAEOR.OH 

AECW?  TOAPIR.O. ROUNO, 

ARIP.3.SF.0L. 

TOACOR. 1 , TOMLR. 2H.SETCS 
AECW  MLP.OH.MRR.IH, 

TOACIR.O. 

ALlP , 3C , ACOR. 1C.APOR.2TS.SF.OH, 


AL=A*DL 

2AL 

F4a-D-Fl»lC 

AM 

3AL*F3*AL. ROUND  AM 

0HaSH.2AL*4AM 

F4 

WRITE  OUTPUT  0 
BRANCH  TO  0101  ON  Fl 
1 AH*F 1 *AH 
AL 

F2*F 1 *D*OC 
OH 

3AL. ROUND  OH 
PASS  AL 
F  2 

AH*A*OH 
1  AH 

F3a-D-F2* 1 C 


50 

' 
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* 

‘ 

* 

» 
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TOMRP.OL 

AL 

AECW? 

RDWRsALlR.  1  .AR.3.4AR, 

READ  SHSL 

MLR. 2L.MPR. OH, ROUND, 

4AHsf4«AH, ROUND 

AL 

ALIR.1L,ACIR.0,ARIR.0,SF. 1H, 

0L=SL*1 AH.3AL 

T0AC0R.2 ♦  T  OMLW. 1L . 

E3 

BN7.P 

BRANCH  ON  CT  TO 

MOPLP 

M0PL1 

AECW 

MLR.?H,MRR.OL, 

?AL*E2*AL 

TOApJR.2. 

AH 

ALl9.3,AC0H,l» aror , ots.se. oh. 

E 1 =E2»D*0C 

toaeor.ol 

OL 

AECW2 

ROW AsMlM. 0 , AL IR.0.AR.2.AAR, 

READ  DHOL 

TOARlR.l, ROUNO , 

„ AH, ROUND  OL 

API  9. 2, SE. GL  , 

PASS  AH 

TOACOR.O.TOMLR.IH.SETCS, 

El 

flCR.q 

BRANCH  TO  0102  ON  E? 

AECW 

SE.OH 

NOR 

AECW  2 

WRM=aEOH.O . AR.4.AAR, 

WRITE  LAST  OUTPUT 

BNZ.C 

BRANCH  TO  STOP 

AECW 

SE.GH 

NOP 

AECW2 

STOP 

STOP 

• 

IN  ro 

OUT  TflANI tions,  u  POSITIVE 

PIOOl 

AECW 

MLR.iH.MRR.OL, 

1 AL»E 1 *AL 

T0APIR.2, 

AH 

ALIR.1H.AR19. 3, SE.1H 

OT=SM.AL 

AECW2 

ROWAsMLR.O,ALIR.O,AR.2,AAR, 

READ  OHOL 

ARIr.2, SE. GL , 

PASS  AH 

toalor.o 

OT 

AECW 

MLR.OU.MRR.IH, 

AL=4*DL 

toapir.o. 

ial 

ACOP.2.SE.0H, 

E4=E  3 

T0MRR. On 

AH 

AECW2 

MLR . ?L , MRR . OL . HOUND , 

3ALsE4*AL, ROUND 

AH 

ALOR.O,ACIR.OC,aRIR.O,SE.1H, 

OH=OT - 1 AH» 1 AL 

TOACOR. 3. T0MLR.2L, 

y  a 

BR.3 

BRANCH  to  M0PL2 

AECW 

MLR. ?H, MRP. OH, 

iah*e?*ah 

T0ARIR.3, 

AL 

ALIR.3.AC0R.1 .SE.OH, 

E2*E?.D 

TOApnR.OH 

OH 

AECW2 

TOABIR.O, 

3AL 

ARIR.3.SE .OL, ROUND. 

PASS  AL, ROUND  OH 

T0AC0R,lfT0ML«.2H 

E2 

Pi  00? 

AECW 

mlr.ol.mrr.ih. 

AL=4*0L 

TOACIH.l, 

?al 

TOMRP.OH 

AH 

AECW? 

ACOP. 2, SE.OH, ROUNO 

F4=E3,  ROUND  AH 

AECW 

MLR. 3L, MRP, OH, 

4AH=f4*AH 

TOAP 1 R , 3 , TO  ACOH , 3 , TOMLR . 2L , 

AL.  E 4 

ACOp.O, SE.OH 

E2sEl 

AECW2 

ARIP.3.SE.0L, 

PASS  AL 

T0AC0R.1» TOMLR. 2H, 

E  2 

SR.7 

BRANCH  TO  M0PL1 

AECW 

MLR, OH, MRP. 1H, 

AH=A*Dh 

T0AP1R.1. 

4AH 

ALIR.3C.AC0R. 1C. AhOR.I , SE.OH. 

E  3=-0-E? .ONE 

TOMRP.OL 

AL 

NSWC  TR  81-313 


P0I01 


P0I02 


MIOOI 


AEC*2 

ROWR=ALlR.l ,Ak.3,AAR, 

PFAU  SMSL 

•73 

MLR.2L.MRP.OH,hOUNO, 

4AH=F  4»AH, ROUND 

AL 

• 

ALlP.lL.ACIR.lC,ARlR.lC,SF.lh, 

OL  =  SL-2AL-<»AH 

• 

TO AC OR. 2*  TOMLR. 1L , 

►  3 

• 

dNZ.O 

BRANCH  ON  CT  TO 

MOPLP 

OUT  TO 

IN  TRANSITIONS.  D  POSITIVE 

AECW 

MLR.0H.MRR.1H, 

ah=a«oh 

*74 

TOACIR.O. 

1  AH 

• 

TOMRR.OL 

AL 

AECW2 

ACOR. 3.SE , OH, ROUND 

F3=F4, ROUND  AL 

75 

AECW 

MLR . 1 L ,MRR . OL  . 

3AL=F3*AL 

*76 

TOARIR.?, 

AH 

• 

TOACOR.2.TOMLR.1L, 

F  3 

• 

ACOR.l.SE.OH 

F1=F2 

AECW2 

RDWA=MLR.0,ALIR.0,AR.2,AAR, 

READ  DHOL 

♦77 

ARIP.2.SE.0L, 

PASS  AH 

« 

TOACOR.O, TOMLR. 1H, 

FI 

• 

HR.  #> 

BRANCH  TO  MIPL1 

AECW 

MLR.OL.MRR.1H, 

AL=A*DL 

*78 

TOAOJR.O. 

3AL 

• 

ALlR.3C.ACOR.OC.AROR. l.SE.OH, 

F  4  =  -D— F 1 .ONE 

• 

TOMRR.oH 

AH 

AECW? 

ROWRsAL IR.1.AR.3.AAR, 

READ  SHSL 

*79 

MLR. 1L.MRR.0L, ROUND, 

3AL=F3*AL. ROUND 

AH 

* 

AL IP. 1L, ACIR.OC.ARIP.OC.SE. 1H, 

OL=SL- 1 4H-3AL 

# 

TO ACOR. 3. T0MLR.2L  . 

F  4 

• 

BNZ.8 

bRANCH  ON  CT  TO 

MIPLP 

AECW 

MLR.2H.MPR.0H, 

2AH=F2*AH 

*80 

T0ARIR.3, 

AL 

• 

ALlR. 1H.ARIR.2.SE ,1H 

OT=SH»ah 

AECW2 

TOALOR.O, 

OT 

*81 

APIR.3.SF.0L 

PASS  AL 

AECW 

MLR. OH, MRP. 1H, 

Ah=A«Dh 

*8? 

TOAPIR.l. 

2AH 

• 

TOMRR.OL, 

AL 

• 

1C0P. 3, SF . OH 

F3=F4 

AECw? 

MLR. 1L.MRR. OH, ROUND. 

4AH=F3*AH. ROUND 

AL 

*83 

ALOR.O.ACIP. 1C.ARIP. 1 .SF.1H, 

0H=0T-24^«2AH 

• 

TO A COR. 2. TOMLR. 1L, 

F  3 

• 

HR.? 

bRANCH  TO  MIPL2 

AECW 

MLR.lH.MRR.OL. 

?J-L-F1*AL 

*84 

T0APIR.2. 

AH 

• 

ALIR.3.ACOR.O.SF.OH. 

E l *F 1 *D 

• 

TOAEOR.OH 

Oh 

AECW2 

ROW AsMLR . 0 , AL 1 R , 0 . AR . 2 • A AR , 

READ  DHDL 

•85 

ARIR.2.SF. OL .ROUND , 

ROUND  OH.  PASS  i 

AH 

* 

TOAPIR. 1 .TOACOR.O. TOMLR. 1H 

El 

IN  TO 

OUT  TRANSITIONS,  0  NEGATIVE 

AECW 

ALIP.IH.ARIR.OC.SF.IH 

Oh=SH-3AL 

86 

AECW2 

TOAEOR.OH, 

Oh 

*87 

ALlP.3C.AC0R. 1C.AR0R.1 .SF.OH 

E4=-Q-F? ♦  1 

AECw 

ALlR, 3, ACOR.l.SE.OH. 

F?sD.F  2 ‘ 

*88 

ml».?h,mrr.oh, 

1 AH=E ?*AH 

• 

TOMLR. 2L, ROUND 

Fa, ROUND  OH 

A-14 
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MI  002 


M0I01 


M0I02 


AECW? 

ALIP.3.AC0P.1.SF .OH, 

TOACOH. 1  * 

BH.7 

F?=D«F? 

F? 

branch  to  mqrl i 

*89 

• 

AECW 

mlr.oh.mrh,.  ih, 

ALlP.3C,ACOH.lC,APOP.i,SF.OH, 

TOACIH.O. 

T0MLR.2H 

AH=A*OH 

F3=-D-F2»ONE 

1  AH 

F  2 

*90 

• 

• 

AECW? 

RUWRsALlH. 1 , AH. 3, AAR, 

MLR. ?L.MHW. OH, 
ALIR.IL.ACIH.O.SE.IH, 

TOACOH. 2. T0MLK.1L, 

BNZ.9 

HEAD  SHSL 

AAHsF  4*AH 

OL=SL* 1  AH 

F  3 

BRANCH  ON  CT  TO 

MOPLP 

*91 

• 

• 

• 

AECW 

SF.AH 

NOP 

9? 

AECW? 

rdwr=mlk.o,alim.o,ar.?.sar 

BACKUP  OHDL  AR.2 

93 

AECW 

ACOR.O.SE.OH 

F?=F  1 

94 

AECW2 

TOACOR. 1 .T0MLK.2H, 

HR. 7 

F  2 

BRANCH  TO  M0PL1 

*9S 

AECW 

MLR.0H.MRR.1H, 

ALlP.3C.ACOR.OC.APGR. 1,SF. OH 

ah=A*Dh 

F3*-D-F1*0NE 

*96 

AECW2 

OUT  TO 

P0WR*ALIR.1,AP.3,AAR, 

MLR. 1L, MRP. OH, 
ALlR.lL.ARIR.lC.SE.lH, 

T0AC0R.2, TOMLR . 1L  » 

8NZ.9 

IN  TRANSITIONS.  D  NEGATIVE 

WEAO  ShSL 

4  AH=F 3* AH 
0L=SL-4AH 

F3 

BRANCH  ON  CT  TO 

MOPLP 

*97 

• 

• 

• 

AECW 

SE.OH 

NOOP 

98 

AECW2 

ARIR, 3C ,SF . OL 

[)LT  =  -Al 

99 

AECW 

T0MLR.3H 

dlt 

100 

AECW2 

ACOP, l.SF.OH, 

HR. ft 

F1*F2 

BRANCH  TO  MIPL1 

•101 

AECW 

MLR.3M,MPR. 1L, 

ALIP.3C.AC0R. 3, SF. OH, 

TOACOH. 0,TOMLP.1H 

AL=-1*ULT 

F  4-F  4-D 

FI 

*102 

• 

AECW? 

MLH.PL.MRR.OL. 

ALIR.1L.AC1R. 1,APIR.0C,SF . lh, 
TOACOR. 3. TOMLW.2L, 

ROWPs  AL I H . 1 , AR . 3 . A AR , 

BNZ.8 

3AL=F4*AL 

0L*SL»?AL-3AL 

F  4 

READ  SHSL 

BRANCH  ON  COUNT 

MIPLP 

*ir«» 

• 

• 

• 

AECW 

T0ARIR.3. 

alir.ih,acir.o,arir.ic,sf.ih 

AL 

OHsSH*! AH-4AH 

*104 

AECW2 

ARIR. 3C ,SF ,0L . 
toafor.oh 

OLTa-AL 

OH 

•10S 

AECw 

MLR.lH.MkP.OU, 

ALIR.3.ACOR.O.SF .OH, 

TOMLR. 3H, ROUND 

1 AL*F 1  * AL 

F 1 =F 1 *0 

DLT, ROUND  OH 

*10ft 

• 

AECW? 

TOACOR. O.TOMLK.IH, 

SR.ft 

FI 

BRANCH  TO  MIPLl 

•  107 

AECW 

MLR.3H,MRR.1L. 

TOAPI R, 0  » 

ALIR.3C.AC0R. 3, SF. OH 

AL=DLT*-1 

1AL 

F  4  =F  4—0 

*108 

• 

AECW? 

P0WBaALIR.l.AH.3,AAR, 

MLR.2L.MWR.0L, 

ALIR.IL.ACIH.O.AWIR.O.SF.IH, 
TOACOR. 3.TOMLR.2L, 

BNZ.A 

READ  SHSL 
3AL*F4*AL 

OL  =  SL ♦ 1  AH, I AL 

F  4 

BRANCH  ON  CT  TO 

MIPLP 

•109 

• 

• 

• 

A-15 
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RSSSSSSS 
SSSSSSSSSS 
SS  S 

rs 

sssssssss 

sssssssss 

SS 

s  SS 

ssssssssss 

SSSSbSSS 


SSSSSSSS 
SSSSSSSSSS 
SS  S 

SS 
SS 
SS 
SS 

SS  s 

SSSSSSSSSS 
SSSSSSSS 


SS 

SS 

SS 

SS 

SS 

SS 

SS 

SS 

SSSSSSSSSS 

SSSSSSSSSS 


SSSSSSSS 
SSSSSSSSSS 
SS  SS 

SS  SS 

SS  SS 

SSSSSSSSSS 
SSSSSSSSSS 
SS  SS 

SS  SS 

SS  SS 


SSSSSSSSSS 

SSSSSSSSS 

SS 

SS 

sss 

sss 

SS 

S  SS 

SSSSSSSSS 
SSSSSSS 


•COMOECK 

SCLA3 

SCL  A  3 

AECW 

Sf.OH 

MOP 

0 

AECW? 

pOWtJ  =  MRR.O.AR.O 

READ  A. fl 

1 

AECW 

SE.On 

MOP 

2 

AECW2 

pOWrtsMpR.l .AP.O.AAR 

RE  AO  c 

3 

AfcCW 

S^.OH 

NOP 

4 

AECW? 

rDwH=MLR.O.AR.  1 

BEAD  XI, XT! 

s 

AECw 

‘tF'.Orl 

NOP 

8 

AECW2 

ROWM=MLW.l«AR.2 

READ  Yl.Y? 

7 

AECW 

MLP.OH.MRk.OH 

AX1=A**1 

8 

AECw? 

pOWb=MLR.2,AR.3« 

Rt  AO  21, Z? 

•9 

AECw 

MLH.OL.MRR.OH 

MLR. IH.MRk.OL.rOAL/P.fl 

AX?=A*X? 

BYi=8*ri.  oes  ax i 

A 

AECw? 

RDWH=ALlR.2.AR.4,SAR, 

DUMMY  read 

SC3LP 

AFCW 

mlr.il.mrr.ol.toalip.i 

MLR.2H.MRW. 1H.TOACIR.O 

fc!Y2=W*Y2.  OES  AX? 

CZ 1 =C*Z1 ,  OES  BY  1 

C 

AECW2 

pDWB=MLR.O.AR. l.AAR, 

READ  X3.X4 

*D 

AECW 

MLR.2L.MRR.1H.T0ACIR.1 

toarir.o 

CZ?=C*Z2.  DES  8Y2 

OES  CZ1 

E 

AECW? 

PDWBsMLR. 1 .AH. 2, AAR. 

READ  Y3.Y4 

*F 

AECW 

TOAWIR.l.SE.OH, 

alIR.o.aciw.o.arir.o 

MLR.OH.MRR.OH.TOAEOR.OH, 

OES  CZ? 

W 1  =AX 1 »8Y 1 *CZ  1 

A  x  3« A  *  X3 ,  OES  W1 

• 

*10 

AECW? 

ALlR.l.ACIR.l.ARlR.l.SE.OH 

powh=mlR.?.ar.3,aar,rouno. 

w?=Ax?.8Y?*CZ2 

READ  L 3, ZA, ROUND  W1 

*11 

AECW 

MLR.OL.MRR.OH, 

SF  . OH . T OAtOR .0L.8N7.9 
mlR.IH.MRR.OL.TOALIP.O. 

AX4»A*X4 

DES  W?,  LOOP 

HY3«H*Y3 »  OES  AX3 

• 

*12 

AECw? 

POUND 

wRMsAEOR.O.AR.A, AAR, 

ROUND  W2 

WRITE  W1.W2 

*13 

AECW 

MLR.lL.MRH.OL.TOALlR.l 
$E  .  OH 

bY4sH*Y4,  OES  AX4 

MOP 

14 

AECW? 

STOP 

STOP 

15 

A-16 
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•COMDECK 

f??S3 


F??LP 


FFFFFFFFFF 
FFFFF  F FFf F 
FF 
FF 

FFFFFF 

FFF FFF 

FF 

FF 

FF 

FF 


F  FF FFFFF 
FFFFFFFFFF 
F  FF 

FF 
FF 
FF 
FF 
FF 

FFFFFFFFFF 

FFFFFFFFFF 


FF  FFFFFF 
FFFFFFFFFF 
F  FF 

FF 
FF 
FF 
FF 
FF 

FFFFFFFFFF 

FFFFFFFFFF 


ffffffff 

FFFFFFFFFF 
FF  F 

FF 

FFFFFFFFF 
FFF FFFFFF 
FF 

F  FF 

FFFFFFFFFF 
FFFFFFFF 


FFFFFFFFFF 

FFFFFFFFF 

FF 

FF 

FFF 

FFF 

FF 

F  FF 

FFFFFFFFF 
FFFFFFF 


F22S3 

AECW  SF.nH 

AECW?  ROWRsMRR.?,  AH. 5 

AECW  SF.nH 

AECW?  RDW0sMWW.3,AR.b,AAW 

AECW  SF.OH 

AECW2  RDWP=MLR.?«AR.l 

AECW  SF.OH 

AECW?  ROWflsMLW • 3.AW.1.AAR 

AECW  MLR.?H.MWW.?L 
AECW2  ROWflsMWH.O. AR. 1 ,AAW , 
MLW.?L.MWW.3L 

AECW  MLR.3H,MRR.?H.T0ACIR.? 

AECW?  ROWPsMLW.O.AR.?. 

MLR.1L.MRR.3H.TOACIR.3 
AECW  TOAptH.? 

AECW?  RCWPsMLW  «1«AR.3« 

TOApIH.3 

AECW  HLR.nH.MWR.OH 

AECW?  ROWRsACIR. 1 .AH.StAAR 

AECW  MLR.IH.MHH.OL.TOALIR.? 

AFCW?  ROWPs AL IW.O.AP.O* 

ALIP.?.ACIR.2,ARIR.?,SF.0C 
AECW  MLR.3M.MHR.2L.T0ALIH.3* 
ALIR.1L. SF.OH, 
T0AC0R.1.T0MRW.2H 
AECW?  WOWPsACIR.OtAR.4.SAR, 
MLR.1L.MWR.3L, 

TOAROR.O. WOUND 

AECw  mlr.ol,hrw.oh,tuawir.?, 

ALlR^ACIR.S.ARlR.S.SF.OL, 

WOUND 

AECW?  ROWHsMLW.P.AR.2, AAW, 

MLR.1L«MKR,OL,TOARIW.3« 
ACOP.l .SF .OH, 

T OApOW • 1 «  T OMWW , 3H 
AECW  MLR. ?H. MRR. 2H, TOAL IR.?* 

ALIP.OH. ACOW. 1 .APOR. 1 .SF.OH. 
TOAEOR. 1H, WOUND 
AECW?  R0WB*MLW.1.A».3.AAW, 

MLR . ?L»MRR.3H,TOALIR.3» 
APOP.l.SF ,0H, 

TOACflW.O,  ROUND 

AECW  MLR .1H.MRW , ?H , 1 OAC 1 W .? « 

ACOP.O.AWOP.OC.SF . 1H, 


NOP 

0 

READ  FH.FL 

1 

NOW 

? 

Rt AD  Gh.GL 

3 

NOP 

u 

READ  A ,C 

5 

NOP 

b 

READ  H.D 

7 

afl=a*fl 

H 

WE AD  U«V 

*9 

CGL*C*GL 

HFM=H*FH,  DES  afl 

A 

WEAO  XH.XL 

*8 

DGH=D*GH.  DES  CGL 

DES  BFh 

C 

WEAO  YH.YL 

*0 

DES  DGH 

UXH*U*XH 

E 

READ  SSH.SSL 

F 

VYH=V«YH.  DES  UXH 

in 

READ  SH.SL 

*11 

FHsUXH* AEL  *BEH 
8FL=w*FL.  DES  VYH 

*1? 

PASS  SSL 

• 

DES  FH 

dummy  read  on  output 

*13 

UGL*D*GL 

• 

DES  SSL,  RNO  Fh 
UXL=U*XL.  DES  SFL 

*14 

gh=vyh*cgl*oGh 

• 

WOUND  SSL 

WE  AD  XH.XL 

•IS 

V YLS V  * YL »  DES  DGL 

• 

fho»fh 

• 

DES  GH 

AFH=A*FH.  OES  UXL 

*lb 

SSH=SH*FH*GH 

• 

DES  FHO.  RNO  GH 

WE AD  YH.YL 

*17 

CGH=C*GH,  DES  VYL 

• 

OHO=GH 

• 

DES  SSH,  RNO  FHO 
c«Fh=H*Fh,  DES  AFh 

•IB 

ShO=SSH ( -SSL ) 

* 

A-17 
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T0AEOR.2H, POUND 

DES  OHO,  RND  SSH 

AECW? 

MLR.1L,MRP.3K«TOACIR.3« 

DGY=D*GH,  DES  C&H 

*19 

ALIp.2.ACIR.?,AR1R.2*SF.0L, 

FL=UXL«AFN*HFL 

• 

TOAFoR. OH, ROUND 

DES  SHO,  RND  OHO 

AECW 

MLR.nH.MRR.flH.TOARIR.?, 

UXM=U*XH,  DES  BFH 

*1  A 

ALlR.3,ACIP.3,ARIR.3,SF.0L, 

GL=VYL*CGH,DdC 

• 

T 0 AroP . ? , T OMRR , 2L  *  ROUND 

DES  Ft,  RND  SHO 

AECW? 

MLR.1H.MRR.OL.TOAR1R.3. 

VYH=V*YH,  DES  D&H 

*18 

ACOp.2,$E .OH, 

FLO=FL 

• 

TOArOR . ? .TOMRR . 3L , ROUND 

DES  GL,  RND  FL 

AECW 

MLR.?H,MHR.?L,TOALIR,2, 

AFL=A*FL,  DES  UXH 

•1C 

ALIP.OL,ACOR.2»AROR.?,SF.OH, 

SSL  =  SL  »FL ♦GL 

• 

TOAf OH. 1L, ROUND 

DES  Flo,  RND  &L 

AECW? 

ROWflsALlR.O,Ak.O,AAR, 

READ  SH.SL 

*1  D 

MLR .  ?L  ,MRR .  3L  .  TO AL I R .  3 . 

LGL=C*GL.  OES  VYM 

• 

APOp.?.SF.OH, 

GLO=GL 

• 

TOApOR.O , ROUNO 

DES  SSL,  RND  FLO 

AECW 

MLR.3H.mRR.?L,T0ACIR.2. 

HFL=H*FL.  OES  AFL 

•IE 

aror.o,acor.oc,sf. iu. 

SLO*SSL (-SSH) 

• 

TOAfoR. 2L. ROUND 

DES  GLO,  RND  SSL 

AEC  W2 

rowr*mlr.o,ar.?,aar. 

READ  XH.XL 

•IF 

MLR.3L,MRR.3L,rOACIR.3, 

DGL=D*GL,  OES  CGL 

• 

ALIr.?.aCIR.2,ARIR.2,SF,0L. 

fh=uxh.afl»rfm 

• 

TOAFOR.OL  » ROUND* BN Z. 8 

OES  SLO.  BNO  &LO,  HR 

AECw 

MLR  . OL ,MRB.  oh.toarir.?. 

UXL=U*XL.  OES  HFL 

*20 

ALlR.3.ACIR.3,ARlR.3.SF.0L, 

GHaVYM«CGL*DGH 

• 

TOACOR . 1 , TOMRR ,2H , ROUNO 

DES  EH,  RNO  SLO 

AFCW2 

wrm=aeor.o,ar.<.,aar. 

WRITE  SHO, SLO 

*21 

mlr.il.mrr.ol.toarir.3. 

VYL=V*YL.  DES  DGL 

• 

acop.i, sf.oh. 

FHO=FH 

• 

TOABOR . 1 , TOMRR . 3H , ROUND 

DES  GM,  RND  Fh 

AECW 

AHOP.O.SF.OM 

PASS  SSL 

2? 

AECW2 

TOAEOR.3L 

DES  SSL 

23 

AECW 

SF.OH 

NOP 

?4 

AECW2 

WRM=aEUR.3,AR,5 

WRITE  SSO 

?s 

AECW 

SF  ,<ih 

NOP 

2* 

AECw? 

WRM=AEUR . ? . AR . 5 » S AR 

WRITE  GO 

27 

AECW 

SF.oh 

NOP 

28 

AECW? 

WRM=aEO«.l,AR.5.SAR 

WRITE  FO 

29 

AECW 

SF.qh 

NOP 

2A 

AECW? 

STOP 

STOP 

28 
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DOOOOODOO  OODDOOOOOD 

DO 

DO  OOOODOOD 

DD 

00 

DODuOOOUOO  OODDOOOOOD 

DODO  DODO  OUUDDDDoD 

ODD 

OD 

DO  00  DO 

DO 

DODO  00  OD  00 

00  0 

DO 

DO  DO  00 

00 

00  DD  00  00 

DD  00 

00 

00  DO  ODODOO 

OD 

OD  DO  DO  00 

00  DO 

00 

DO  DO  OODDOD 

OD 

DO  DD  DO 

OD  DO 

00 

DO  00  DO 

DO 

DD  DD  00 

DO  00 

00 

on  on  oo 

OD 

DD  OD  00 

DO  0 

DO 

OOODOOODOO  DDOODDOODU 

DO 

DD  UDOOUOOCOO 

00 

ODD 

DOOUUOOOU  DOODDOODDO 

00 

DU  ODODDDOD 

OD 

00 

•COMOECK  DEMON 

DEMON  AECW 

SF  ,0H 

NOP 

0 

AECW? 

R0WRsALlR.2.AH.(j 

READ  Ll.Tl 

1 

AECw 

SF.OH 

NOP 

2 

AECW2 

PDWR  =  MLW.n,  AR.fi.AAR 

READ  K 1 , A 1 

3 

AECW 

ALIR.2L.SF.1L 

PASS  Tl 

4 

AECW? 

ROWRsMWR, 1.AR.6.AAR, TOACOR . 0 

READ  H.C.  OES  Tl 

5 

AECW 

SF.OH 

NOP 

6 

AECW2 

ROWP*AL IR.3.AR.6.AAR 

READ  L2.T2 

7 

AECW 

sf.oh 

NOP 

8 

AECW? 

R0WR=MRR.2.AR.4, 

READ  a  11. *12 

•9 

SCMR.OH=AR.ft 

cos  ism  aeo 

AECW 

ALIR.3L.SF.  1L 

PASS  T 2 

A 

AECW? 

R0WB*mLR.2«  AR. 7 .TOACOR. 1, 

READ  Y11.Y1?.  DES  T? 

*B 

SCMR.0H=AR.2 

COS  (Sll>  AE1 

AECW 

sf.oh 

NOP 

C 

AECW? 

R0WR=MLR.l,AR.6.AAR 

READ  K?,A? 

D 

AECW 

SF.OH 

NOP 

E 

AECW2 

RDWfl=MHH.3,AR.4,AAH, 

RE AO  A21.X22 

*F 

SCMP.OHsAR.l, 

COS  (S?l>  AEO 

• 

MLR.f)L,MRR.?H 

A* 1 1 sA 1**11 

AECW 

MLR.?H,MRR. 1L 

cyh=c*yi  i 

10 

AECW? 

ROWRaMLR. 3 , AR. 7 . AAR. 

head  Y21.Y2? 

*1  1 

SCMP.OHsAR.3. 

COS  (S?l>  AEl 

• 

MLR. ?L.MRR. 1H.TOARIR.O 

HY 1 ?®R*Y 1 ? .  CtS  AX11 

AECW 

mlr.oh.mrr.oh.toacir.o 

KS11*K1*S11.  DES  CY 1  1 

12 

AECW2 

TOALIR.O 

DES  PY1? 

13 

AECW 

T0APIR.2 

OES  KS11 

14 

AECW? 

ROWPsALlR. 1 .AR.5. 

Pt AO  TOl.TG? 

•IS 

SCMP.OHsAR.OlNC. 

COS  (SI?)  AEO 

* 

MLR.3H.MPP. 1L 

CY?1»C*V21 

AECW 

MLR. JH.MRP.OH. 

KS21=K2*S2l 

*16 

AlIR.?H,ACOH.OC.ARIP.?,5F.1H 

P11=L1-T2*KS1 1 

AECW? 

ROWflsAC 1R.3.AH.5.SAR, 

Dummy  READ  ON  OUTPUT 

*17 

SCMP,OHsAR.?INC, 

COS  (S12)  AEl 

• 

MLR. 3L. MRP. IH.TOACIR.1. 

tsY??*B*Y??.  OES  CY?1 

* 

ALlP.O.ACIR.O.ARlR.O.SF.OL, 

Yl 1=RY1?*CY1 1*AX1 1 

• 

TOAPOR.O.SETCS 

OES  Pll 

AECW 

MLR.1L.MHR.3H.T0AR1P.?, 

AX21sA2*x?l,  OES  KS21 

*  1  8 

ACOP . 0 . AHOR. OTb.SF .1L, 

0M1 1  =  T 1  * (PI 1)C 

• 

TOALOP . 0 . TOMLR. 2H 

DtS  Yll 

AECW? 

SCMP.OHsAR. 1  INC. 

COS  ( S?2>  At' 0 

•19 

MLR. ?L. MRP. 1L. TOALIR.O, 

CY 1 2=C*Y1 2,  OES  9Y2? 

• 

TOMBP.OL.EOUNO 

OtS  DM1 1 ,  PNO  Yl  1 

AECW 

MLR. 3L ,MkP. IL.TOAnIR.}, 

CY??*C*Y22«  DES  AX?  1 

*  1  A 
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ALiP.3H,ACOR.lC,ARIR.2,SF.lh 
AECW2  SCMP.OHsAP.UnC, 

MLR.?H,MRR.OL,TOACIR.O, 
ALlP.O.AClP.l.ARiW.l .SF.OL, 
TOApow.O.SETCS 

AECW  MLK.OH.MRP.OH.TOACIR.l, 

A COP, 1 ,  AROB  .0TS.SF.1L, 

TOALOR. 1 . TOMLR.3H 
AECW?  Mt-R.?H.MRR.lHtTOAClR.2» 
ALOp.O.SF.OH, 

TOMpp . OL .POUND 

AECW  mup.oL,HPR.2L.TOARIR.2, 

TOAEOR.OH 

AECW2  RDwPsMRR.2. AP.4,AAR, 
SCMp.0HsAR.0lNC, 
MLR.3H,MKR.OL,TOALIR.O. 

ALlR.2H.ACOR.0L,ARIP.2,SF. 1H 

AECW  MLR. 1H,MRP.OH,TOAHIP.O, 

ALOP.l  .SF.OH,TOAROR.O,StTCS 
AECW?  SCMP.0HsAR.2lNC, 

MLR . 3H  «  MttR . 1 H , TO AR I P . 3 . 
ALlR.0.AClR.0,ARlR.0,SF.0L. 
TOAEOR- 1H 

AECW  MLR. iL,Mnr.3L,TOAR1R.?, 

acop.o, arop.ots.sf.il, 

TOALOR.2.TOMLH.2L 
AECW2  SCMP.OHaAR.llNC, 

MLR.pH.MRR. 1L.T0AL IP. 0. 
ALIR.3H.AC0R.1C,ARIP.?,SF.1H, 
TOMoo.OL .POUND 

AECm  MLR.Ifl.MHW.lL.TOARIR.l, 

ALIP.1H.ACIR.2.ARIR.3.SF.2H, 
TOApnR.O.SETCS 
AECW2  SCMP.0H=AP.3lNC, 

MLR.pL ,MRR .OL.TOACIR.O. 

ALIR.O.ACIR.I.APIP.I.SF.OL. 

T0AF0R.2H 

AECW  MLR.flH.MRR.OH.TOACIR.l. 

ACOR. I .AROR.OTS.SF.1L, 

TOAl  OR . 3 , TOMLH . 3L . ROUND 
AECW?  MLR.PL.MHR.IH.TOACIR.P, 
ALOP.2.SF.OH, 

TOMpp.OL, ROUND 

AECW  MLR.0L.MRR.2H.TOARIR.2, 

TOAEOR.OL 

AECW?  ROWRsMPR. 3.AH.A.AAR, 
SCMp.OHsAR.OlNC, 
MLR.TL.MRR.OL.TOAlIR.O, 
ALIR.2H,ACOR.oC.ARIR.2,SF.IH 
AECW  MLR . ] H.MWR .OH.TOARlR.O. 

ALOP . 3 . SF . OH , TO AROH . 0 . SE  TC S 
AECW?  5CMR.0HSAR.2INC, 

MLR. OL.MRR.lH.TOARl P.3, 
ALIR.O.ACIP.O.AHIP.O.SF.OL, 
TOAEOR. 1L 

AECW  MLR.lL,MHR.3H,T0APJR.2, 

A COP, 0 , APOR, OT  S, SF . 1L  » 

TOAl  OR.0.TOMLM.2H 


P21=L2-T2*KS21 
LOS  ( 52? •  AE1  *1« 

01  1  s|)M  1  1  *Y1  1  ,DES  CY  1 2  * 
Y21sBY??,CY21*A*21  * 

DFS  P?1 

ES1?=*'  1*S1  ?,  OES  CY22  *1C 
DM?] =T2* (P21 >C  • 

OtS  Y?l 

BY  1 1 =H* Y l l .  DES  Dll  *10 
PAS  Y 1 1  * 

DES  Dm? i ,  UNO  Y21 
AX  1 ?=A 1 »X 12 »  DES  «S12  *  IF 
DES  Y 1 1 

READ  XU. XI?  *1F 

COS  (Sll)  AEO  • 

DPI sDM?l • y2 l ,DES  BYll  • 

P 1 2=L 1 - T 1 »KS12 
KS??»K?*S?2,  DES  AX1 2  *20 
PASS  Y? 1 .  DES  PI 2 
COS  IS11)  AE1  *21 

P Y? 1 =8* Y2 1 ,  DES  021  * 

Y12=BY1 1 »CY12*AX1?  * 

DFS  Y? 1 

AX??=A2*X22,  DES  KS22  *22 
DMl?sTl. (P121C  * 

DES  Y12 

COS  IS?1>  AEO  *23 

C Y 1 1 =C*Y 11.  DES  BY21  • 
P2?sL2-T?.KS22  * 

OtS  DM]P.  RNO  Y12 
CY?1*C*Y?1.  DES  AX??  *?4 
TDlsTGl.Dll.021  * 

OFS  PP? 

CUS  <  S? 1 )  AE1  *2S 

D12sOMi?»Yl?,OES  CYl 1  * 
Y2?=BY?1*CY?2*AX2?  * 

DES  TGI 

K S 1 1  sK  1  *S1  L  .  OES  CY21  *?f> 
Dm??sT2. (P221C  * 

UES  V??,  RNO  TGI 
PY1?=H*Y1?.  DES  D12  *?7 

PASS  Yl?  • 

DES  DM2?,  RNO  Y22 
AXl lsA\*Xl  1  .  OES  KS11  *?S 
OES  Yl? 

READ  X?1 ,X??  *?R 

COS  (Si?)  AEO  * 

0??=0M22«Y?2.DES  BY  1 2  * 
P11=L1-T1*KS11 
«S21*K2*S21.  DES  AX  1 1  *?A 
MASS  Y22,  OtS  PI  1 
COS  ( Si  2)  AE1  *2B 

HY??sM*y?2,  DES  022  • 

Yl 1=BY1?*CY1 1»AX1 1  * 

DES  Y2? 

AX21SA2*X21,  DES  KS21  *?C 
DMllsTl.(PU>C  • 

DES  Yll 
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AECW2 


AECW 


AECW? 


AECM 


AECW? 


AECw 

AECW? 

AECW 

AECW? 

AECW 

AFCW? 


P0WR=ALIR.l,Aw.S.AIRO, 

SCMP.0H=AR.  1  INC, 

MLR.?L,MHR.1L,  tOAUR.O* 

ALlR.3H,ACOR.lC,ARIR.2,SF.lh, 

TOMQP.OL, round 

MLR.U,MRR.1L,TOARIR.i. 

AUp.1L,ACIR.?,ARIR.3,SF.?m. 

TOApOH.O  * St TCS 

SCMP.0H=AR.3INC, 

MLR.?M,HRR.OL,TOACIR.O, 

ALIP.0,AC|R.1,ARIR.1,SF.0L, 

TOAFnR.2L,PNZ.K 

MLR.flH.MKR.OH.TOACIR.l, 

ACOP.1,AROR.OTS,SF.R, 

TO AL  0« . 1 , TOHL w . 3H , ROUND 
WRM*AtOH.2.AR.S,SAP, 
MLR.?H,MKR.1H,T0ACIR,2, 
alop.o, sf.oh, 

TOMPp.OL, ROUND 
SF.oH 

WRM=AEOR.O,AP.7,SAR 

SF.OH 

WRMaAEOH «1,AR.7«AAk 

SF.OH 

STOP 


READ  T61,TG?  *?D 

COS  (S?2)  AEO  • 

CY l 2=C*Y1 2,  DES  9Y2?  * 


H21=L?-T2»KS21  * 

DES  DM11,  RNO  y 1 1 
CY??*C*Y22,  DES  AX21  *2E 
TG?=TG?»0l2*D2?  * 

DES  P?1 

COS  CS2?)  AE1  *?F 

Dll *OM 1 1 ♦ Y 1 1 ,  DES  CY12* 
Y?1=BY22*CY?1*AX21  * 

DES  TG?,  BR  TO  DMNLP 
KS12=Ki.S12,  DES  CY2?  *30 

DM?1*T2, (P21 1C  • 

DES  V?l,  RND  TG? 

WRITE  TG1.TG2  *31 

Ryh*h*Y11.  OES  DU  * 

PASS  YU  • 

DES  DM?l,  RNO  Y21 
NOP  32 

WRITE  STATES  V11.Y12  33 

NOP  34 

WRITE  STATES  Y?1 , Y22  35 

NOP  36 

STOP  37 
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